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Channel  Modeling  and  Threshold  Signal  Processing 
in  Underwater  Acoustics:  An  Analytical  Overview 

1.  Introduction 

The  purpose  of  this  paper  is  twofold:  (1)  to  combine  a  canonical 
characterization  of  random  underwater  acoustic  fields  or  "channels"  with 
the  principle  elements  of  threshold  signal  processing,  in  essence,  signal 
detection  and  extraction;  and  (2)  to  provide,  in  the  process,  a  necessarily 
concise  analytic  description  of  some  of  the  principal  formal  approaches 
required  to  achieve  these  ends.  Accordingly,  our  presentation  is  primarily 
tutorial,  in  keeping  with  the  space  available.  This  overview  seeks  to 
provide  one  possible  albeit  rather  formal  framework  for  handling  the 
increasingly  specialized  problems  of  medium  modeling,  signal  processing, 
and  performance  evaluation,  which  are  now  encountered  in  applications.  We 
shall  lean  heavily  on  some  of  the  author's  recent  work  on  channel  charac¬ 
terization  [l],  [2]  and  space-time  signal  processing  [3]-[5],  as  well  as 
on  recent  related  work  on  the  stati stical -physical  modeling  of  nongaussian 
noise  and  interference  fields  [6]-[ll].  (Tutorial  background  studies  for¬ 
mally  relevant  to  the  above  are  provided  in  [12],  [13]  for  nongaussian 
electromagnetic  interference  (EMI)  environments.)  In  large  part  we  shall 
be  '-oncerned  with  optimal  processing,  considered  as  a  limiting  form  of 
system  structure  and  performance,  to  be  approximated  in  practice  according 
to  the  inevitable  economic  and  operational  constraints  of  the  task  at  hand. 
[For  related.,  suboptimum  systems,  <*•«  note  the  parallel  telecommunication 
examples  descussed  in  [4],  [14]-[16].] 

Thus,  our  general  aim  is  to  provide  a  framework  for  integrating  the 
effects  of  the  channel  (or  medium)  and  transmission  and  reception  of  sig¬ 
nals  through  it,  to  obtain  system  optimization  and  associated  performance 
measures.  The  complex  acoustic  channels  specifically  considered  here  are 
the  essentially  linear  underwater  media  of  typical  ocean  environments, 
which  include  volume  and  interface  (surface  and  bottom)  scatter  mechanisms, 
and  other  inhomogeneities  (e.g.,  gradients,  internal  waves,  etc.).  Ambient 
noise  mechanisms  must  also  be  considered.  These  include  shipping,  bio¬ 
logical  phenomena,  and  various  geophysical  sources,  e.g.,  arctic  ice. 
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seismic  activity,  etc.  For  example,  see  the  recent  work  of  Dwyer  [17  and 
refs,  therein]  on  underice  ambient  noise.  For  the  signal  processing  fol¬ 
lowing  reception  we  are  mainly  concerned  here  with  threshold  operation,  up 
to  the  point  of  signal  decoding  procedures. 

Because  the  channel  exerts  such  a  critical  influence  on  the  signals 
transmitted  through  it,  it  is  particularly  important  to  include  the  princi¬ 
pal  effects  whereby  these  non-ideal,  inhomogeneous  media  degrade  and  modify 
the  original  signal.  A  canonical  characterization  [namely,  one  that  is 
invariant  of  any  particular  medium  properties  (until  the  medium  is  specified)] 
is  needed  in  order  to  provide  the  required  generality  in  compact  and  manage¬ 
able  form,  and  to  allow  appropriate  approximations  in  specific  cases.  Simi¬ 
larly,  it  is  well  understood  that  effective  signal  processing  must  incorporate 
the  relevant  characteristics  of  the  channel,  as  well  as  the  desired  signal, 
optimally  for  limiting,  optimal  performance  and  at  least  adequately,  for 
practical  suboptimum  systems,  which  may  be  more  or  less  close  to  the  ideal, 
optimum  limits. 

A  critical  element  of  signal  processing  is  spatial  sampling,  achieved 
through  the  array  or  aperture  by  which  source  and  receiver  are  coupled  to 
the  medium.  In  addition  to  the  usual  temporal  processing,  spatial  processing 
becomes  important  in  reducing  minimum  detectable  signal  levels  whenever  the 
interfering  noise  fields  are  nonuniform  over  the  array  or  aperture,  as  long 
as  signal  wavefront  levels  are  maintained  uniform  [3]-[5].  This  is  a  further 
reason  for  examining  the  field  structure  of  the  acoustic  channel  (or  any  chan¬ 
nel,  for  that  matter),  and  for  including  the  specific  effects  of  the  coupling 
arrays  in  tlie  overall  processing  program.  A  second  critical  feature  is  the 
practical  fact  that  the  relevant  statistics,  i.e.,  the  various  probability 
distributions  which  describe  the  noise  fields  and,  in  part,  the  desired  • 

signal  fields  as  well,  are  very  often  highly  nongaussian,  characteri Stic  of 
ship,  underice,  and  biological  mechanisms,  among  others  [3],  [5].  A  third  t 

important  concept  is  that  of  threshold  signal  processing,  where  it  is  pos¬ 
sible  to  obtain  canonical  algorithms,  whose  forms  are  independent  of  the 
particular  noise  statistics  and  signal  structures,  and  which  for  suitably 
weak  signals  are  both  locally  and  asymptotically  (i.e.,  as  processing  gain 


2 


is  increased)  optimum  [4],  [13],  [16].  Applying  the  specific  (non)gaussian 
distributions  indicated  by  these  algorithms  then  yields  the  corresponding 
processing  and  performance.  In  practice,  of  course,  optimal  structures 
are  limiting  forms  which  are  only  approximated  practically.  However, 
using  the  theoretical  weak-signal  optimum  with  care  to  avoid  destruction 
of  the  i nformation-bearing  portion  of  the  signal,  provides  an  effective 
guide  to  suboptimum  algorithms  which  are  "close  to"  the  optimum  and  which 
perform  effectively  at  all  signal  levels:  [18],  for  example. 

Because  space  here  limits  the  detail  with  which  we  can  accomplish 
our  dual  task  of  channel  characterization  and  threshold  signal  reception 
for  these  underwater  acoustic  regimes,  our  treatment  is  basically  a  "top- 
down"  approach,  proceeding  from  rather  general  formulations  to  more  specific 
examples  which  allow  us  to  invoke  some  of  the  physics  of  the  problems  in 
question.  [Extensive  references  assist  the  reader  in  his  pursuit  of  ana¬ 
lytical  detail  and  special  solutions.]  One  advantage  of  the  "top-down" 
approach  is  that  it  readily  reveals  the  broad  connections  between  channel 
physics,  effective  signal  processing  procedures,  and  the  general  communi- 
cational  tasks  of  signal  detection  and  estimation  in  noisy  environments , 
focused  here  on  underwater  acoustic  milieux. 

1.1  Organization 

Accordingly,  the  present  paper  is  organized  as  follows: 

Section  2  gives  an  introductory  review  of  the  general  problem,  including 
the  channel  and  the  generalized  communication  processes  involved. 

Section  3,  beginning  Part  I:  Medium  Modeling,  discusses  some  typical 
propagation  equations  for  inhomogeneous  deterministic  and  random 
media  (Langevin  equations)  and  their  formal  solutions.  This  includes 
a  number  of  recent  results  developed  for  dealing  with  (linear) 
random  media,  in  particular  complex  scattering  channels  like  the 
ocean,  which  combines  a  variety  of  interacting  scattering  mech¬ 
anisms.  [An  equivalent  diagrammatic  formulation  is  outlined  in 
Appendix  A.  1.] 

Section  4  presents  an  outline  of  a  second-order  statistical  representa¬ 
tion  of  these  channels,  while 
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Section  5  gives  two  important  illustrative  examples  for  underwater 

acoustics:  I,  weak  volume  scatter,  and  II,  an  "exact"  ocean  surface 
scatter  formalism,  with  a  short  discussion  of  their  physical 
properties  and  structure. 

Section  6  concludes  Part  I  with  a  summary  of  poisson  field  statistics 
which  are  needed  in  the  description  of  both  ambient  and  scattered 
fields. 

Sections  7  and  8,  introducing  Part  II:  Acoustic  Threshold  Signal  Pro¬ 
cessing,  summarize  some  of  the  main  results  for  threshold  signal 
detection,  including  space-time  sampling,  performance  results, 
and  examples  of  optimum  and  suboptimum  detection. 

Section  9  then  provides  a  short  summary  of  signal  estimation  pro¬ 
cedures,  which  are  closely  related  to  the  detection  formalism  in 
the  limiting  optimum  cases.  As  in  detection,  these  are  likewise 
applied  to  threshold  signal  operation,  including  a  specific 
example  of  weak  signal  estimation. 

Section  10  completes  our  analytical  overview  of  signal  processing 
problems  associated  with  the  underwater  acoustic  environment, 
including  supporting  references  to  the  technical  procedures 
presented  in  the  accompanying  References. 

[Supplementing  the  above,  Appendix  A.l  outlines  elements  of  diagram 
methods  used  here  for  the  solutions  of  scattering  problems,  while 
Appendix  A. 2  provides  a  very  concise  summary  of  decision-theoretical 
formalism  needed  in  the  signal  estimation  aspects  of  reception.]  Various 
new  results,  appearing  in  the  cited  references,  are  briefly  summarized  when 
appropriate.  We  stress  once  more  the  isomorphism  between  the  channel 
modeling  and  signal  processing  structures  of  the  (scalar)  forms  arising 
for  telecommunications  in  EMI  environments  and  that  reviewed  here  for 
underwater  acoustic  applications.  Finally,  we  emphasize  that  what  is 
ultimately  needed  for  an  effective  treatment  of  the  general  processing 
problem  is  an  interdisciplinary  approach,  ranging  from  the  appropriate 
physical  models  to  the  methods  of  statistical  communication  theory  (SCT). 

This  will  be  amply  demonstrated  in  the  following  sections. 
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2.  Formulation 

The  general  situation  is  embodied  in  the  operational  relation  [2], 

[19] _ 

tv)  •  •  <2-)) 

where  { u }  is  a  set  of  "messages"  to  be  transmitted  and  { v }  is  the  ensemble 

of  received  messages,  or  decisions,  which  are  consequent  upon  the  set 

{ u } .  TheJ^-p  T^  ,  JAR  are  general  operators  which  describe  respectively 

the  process  of  coupling  to  the  medium  by  the  transmitting  aperture,  the 

effects  of  the  medium  itself  on  the  resulting  injected  space-time  signals, 

and  the  receiving  aperture,  by  which  the  received  field  is  converted  into 

a  (temporal)  wave,  for  further  signal-processing  by  the  receiver  itself. 

The  operators  Tt  ,  T„  are  temporal  processes  only,  representing  the 
i  o  ~‘Ko 

overall  "encoding"  and  "decoding"  of  the  original  and  final  "message" 
sets,  {u},  {v},  in  (2.1).  This  Tt  includes  the  actual  encoding  process 
whereby  "messages"  are  converted  to  signals,  and  then  are  suitably 
modulated  (as  narrow-band  waves,  usually),  to  drive  the  transmitting 
aperture  J[AR.  Conversely ,  includes  the  corresponding  decoding  process, 
as  well  as  any  appropriate  (usually  non-linear)  signal  processing  [20], 
as  shown  in  Fig.  2.1. 

In  our  present  study  we  are  concerned  essentially  with  the  "signal" 
portions  of  (2.1),  namely  the  injected  signal,  Sip(t),  the  received  wave, 
X(t),  following  the  receiver  aperture  JAR,  and  ultimately  here,  with  the 
detected  wave  Y  =  Ioet{X}  (before  any  decoding),  viz.: 

(2.2) 

where  S-n  is  an  injected  (possibly  encoded)  signal  and  where 

(2.3) 

is  the  field  at  some  point  (R_,t)  in  the  medium  in  question,  and  is  a 
source  function.  Here  Gy  and  X(t)  are  explicitly  given  by 


5 


,a' 


1 u  >o—l  T, 


mod 


Encoding 


Trans.  Aperture 


Medium 

Operator 


Rec.  Aperture 


^AT 


S1n^) 


■/v(hT.ar)dvT 


-Generalized  Channel 


kW  T. 


b' 


Oet 


hjH  Tdemod 


of  v> 


''Oetectiorj'1  Oecodlng 
V(t) 


Figure  2.1  Operational  schematic  of  a  generalized  channel:  here,  a 

complex,  ocean  medium  and  subsequent  signal  processing  (TD  ). 


GT  =  GTSin  *  =  /  hT(t-T,tU)Sin(T,5)dT  (l  lySJ, 

(2.4a) 

and 


X(t)  =  Ra  =  /  dVD(R)  /  hD(t-T,t|ReVD)a(R,T)dT  =  J  hR*adVR! 

R  R  ~  —  R  ~  R  ~  RK  K 


(2.4b) 


where  the  source  and  array  operators,  G-p  R  are  respectively t 

CO 

Gt  =  GT(t,JJx)  =  /  hT(t,t||)  (  )Rjt _Tdx  ; 


(2.5) 


R  i  R(t|R.T)  -  /  <)VR  S  h  (T,t|REVR)  (  )R  ,TdT  (■  Tar).  (2.6) 
R  -00  — 


We  remark  that  h^  R  are  respectively  the  filter  weighting  functions, 
now  four-dimensional  (i.e.,  spatio-temporal)  quantities,  of  the  transmitting 
and  receiving  apertures,  which  may  be  time-varying  as  well  as  frequency 
selective.  Their  space-time  fourier  transforms  are  the  (complex)  beam 
patterns 


a 

T,R 


(v,f 


t) 


=  n:\{hT  p 

~t  T,R 


>  *  / 


dTe 


2-n  i  v 


R,V, 


v*Cr 


“  1  CaJT  i 

hj  ,Re 


2rrf  , 

(2.7) 


with  v  a  vector  spatial  frequency  (which  may  itself  be  a  complex  function 
of  frequency,  f),  when  the  medium  is  absorptive. 


tThroughout,  we  denote  operators  by  (  ). 
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2.1  The  Generalized  Channel 

For  the  linear  media  assumed  throughout,  the  field  a  obeys  a  partial 
integro-differential  equation  of  the  form  [2] 

(2.8) 

WQ)  '"(01 

where  L'  '  is  a  linear  (scalar)  partial  differential  operator  with  Lv  ' 
associated  with  the  homogeneous  portion  of  the  medium.  Here  Q  is,  in 
general,  a  (scalar,  linear)  integro-differential  scattering  operator, 
which  describes  the  interaction  of  the  incident  (homogeneous)  field  with 
the  (differential,  or  local)  scattering  (i.e.,  re-radiating  elements)  of 
the  inhomogeneous  portion  of  the  medium.  Boundary  (b.c.)  and  initial 
conditions  (i.c.)  are,  of  course,  a  necessary  part  of  (2.8),  as  indicated. 

Since  we  are  dealing  generally  with  random  media ,  inasmuch  as  Q 
contains  random  as  well  as  deterministic  components,  we  are  concerned  with 
the  appropriate  ensembles,  {a},  (X),  etc.  The  ensemble  of  equations  (2.8) 
governing  propagation  of  the  field  is  now  the  Langevin  equation  ([2];  [20], 
chapter  10).  Since,  formally, 

(a  =  (L(0)  -  Qf1  (-Gt)}  ,  (2.9) 

comparison  with  (2.3)  shows  that 

{T^N)  =  (L(0)  -  Q)'1} .  (2.10) 

This,  and  Eqs.  (2.2),  (2.3),  are  shown  schematically  in  Fig.  2.1. 

For  the  purely  deterministic  special  cases  we  seek  solutions  of 
(2.8)  directly.  However,  for  the  general  situation  of  random  media,  the 
"solutions"  of  the  corresponding  Langevin  equation  (2.9)  are  the  various 
statistics  of  the  field  {a},  e.g.,  the  various  moments  <x>,  etc., 

and  the  moments  of  the  (linear  functionals)  of  the  field,  { X ( t ) > ,  cf. 

(2.4b),  e.g.,  <X>,  <X^X,p>,  etc.,  and  more  comprehensively,  the  various 
probability  densities  (and  distributions)  of  a  and  X.  Thus  we  are  inter¬ 
ested  here  primarily  in  the  means  <a>,  <X>,  the  intensities  ,  <X^>, 
and  the  covariances  (-.Cu^X  -^X^).  [Appendix  A.l  provides  a  short 
introduction  to  various  methods  of  evaluating  <£>,  ,  etc.]  Accordingly, 
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for  the  purposes  of  this  paper,  we  shall  define  the  generalized  channel 
here  as  the  sequence  of  operations  indicated  by  (2.2),  viz.  X(t),  or  in  a 
less  restricted  sense,  by  the  resultant  field  a,  (2.3),  cf.  Fig.  2.1. 
Consequently,  the  desired  statistical  description  of  the  generalized  channel 
is  just  that  of  X  =  Ra,  or  of  a,  namely,  the  various  above-mentioned  moments. 

When  X  is  a  gaussian  process,  these  first-  and  second-order  moments  are 
sufficient  for  subsequent  signal  processing  (T^  in  Fig.  2.1),  to  yield  { v } . 
But,  as  is  often  the  case,  X  here  is  not  gaussian,  so  that  at  least  first- 
order  distributions  of  X  must  be  developed  for  effective  processing.  This 
is  discussed  further  in  Section  6.2  below. 


2.2  The  Role  of  Scattering  Theory 

The  stochastic  character  of  the  medium  or  channel  (as  above),  arises, 
of  course,  from  the  random  spatio-temporal  inhomogeneities  within  and 
bounding  the  medium.  Our  quantitative  description  of  the  medium  now  as  a 
generalized  communication  channel  accordingly  must  be  made  in  terms  of  an 
appropriate  scattering  theory,  particularly  one  which  is  capable  of  in¬ 
cluding  both  the  geometrical  and  statistical  features  of  real-world  situ¬ 
ations  in  manageable  fashion.  We  include  formally  as  a  special  case  of 
scattering  models  that  of  ambient  noise,  where  now  the  secondary  radiating 
or  scattering  sources  are  replaced  by  primary  sources,  cf.  Sec.  6.1. 

Some  recent  new  approaches  [l]  provide  additional  technical  and 
insightful  methods  for  "solving"  the  Langevin  equations  (2.9)  of  this 
general  underwater  acoustic  channel.  These  include: 

I.  The  concept  of  Feedback  Operational  Solutions  (FOS),  cf.  Sec.  3.2 
ff.,  to  represent  the  canonical  Langevin  equation  (2.9); 

II.  The  ability  to  handle  possibly  interacting  scattering  modes 
(the  so-called  "M-form")  (Sec.  2.5  of  [1]);  and 

III.  The  observation  that  boundaries ,  wave  surface  and  ocean  bottoms 
for  example,  can  be  treated  as  distributed  inhomogeneities 

in  an  otherwise  infinite  medium  (Sec.  2.2,  [1]). 

The  first  item  (I)  relates  the  analytic  solutions  and  diagram  methods  to 
a  generalized,  i.e.,  four-dimensional,  feedback  theory  and  formalism; 
the  second  (II)  offers  a  simple  taxonomy  for  accounting  for  any  significant 


scattering  interaction,  e.g.,  surface  with  volume  scatter,  etc.;  and  the 
third  (III)  greatly  simplifies  the  analysis,  at  least  in  the  usual  far- 
field  (or  "radiation-field")  cases,  by  effectively  converting  distributed 
(or  integral)  boundary  conditions  into  local,  or  differential  boundary 
conditions,  i.e.,  plane-wave,  at-a-point  b.c.'s. 

Finally,  we  must  note  that  scattering  theory,  particularly  with 
respect  to  acoustic  and  electromagnetic  propagation,  is  a  venerable  field 
of  study,  with  a  correspondingly  vast  literature.  A  full  citation  of  the 
principal  publications  to  date  is  well  beyond  our  capabilities  here. 

However,  a  useful  and  extensive  summary  of  modern  "classical"  methods  and 
results,  chiefly  in  the  electromagnetic  cases,  has  recently  been  given 
by  Ishimaru  [22],  along  with  extensive  references.  An  earlier  and  more 
mathematical  discussion  of  scattering  problems  has  been  provided  in  the 
interesting  review  paper  by  Frisch  [23].  For  acoustical  scattering,  par¬ 
ticularly  in  the  ocean  and  at  its  interfaces,  the  somewhat  later  review 
articles  by  Fortuin  [25]  and  Horton  [26]  regarding  surface  scatter  are 
especially  noted.  In  the  latter  connection,  see  also  [11],  [43] -[56] ; 
in  particular,  Bass  and  Fuks  [48]  and  [43]-[45], [47] ,  [49],  and  [50],  for 
wave-surface  scattering.  Other  important  recent  studies,  devoted  prin¬ 
cipally  to  volume  scatter,  are  described  by  Flatte  et  al .  [46].  Tatarskii's 
significant  work  [24]  is  also  well  known  and  is  also  especially  appropriate 
to  the  treatment  of  both  acoustic  and  electromagnetic  scattering  in  the 
atmosphere,  and  contains,  as  well,  many  useful  references. 

An  important  departure  from  "classical  continuum  theory"  [25],  [26] 
is  the  so-called  FQM  theory  (after  Faure  [28],  Ol'shevskii  [29],  and 
Middleton  [31].  This  latter  is  based  on  first-order,  i.e.,  independent 
poisson  point-scatter  models.  Various  elements  of  the  author's  more 
recent  approaches  [2],  which  include  both  the  "classical"  and  FOM  theories 
as  particular  components  of  the  general  multiple-scatter  regime,  cf. 
Appendixes  A.l,  A. 2  have  already  been  presented  in  different 
stages  of  evolution  and  development  since  1973,  [32],  [40].  The  essen¬ 
tially  mature  form  is  described  in  [2]. 
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Part  I.  Medium  Modeling 


3.  Inhomogeneous  Stochastic  Media: 

—  — ,i  i  "  'in  — «  -  — ..I— ^  i 

Langevin  Equations  and  Solutions 


Starting  with  the  appropriate  equations  of  state,  mass,  and  momentum 
conservation,  one  can  obtain  by  suitable  approximations  ([2],  [21])  of 
the  resulting  nonlinear  propagation  equations,  1 ineari zed  "wave  equations, 
which  are  special  variants  of  (2.8).  It  is  more  convenient,  however,  to 
proceed  with  the  canonical  form  (2.8)  directly.  Here  is  a  linear 

A 

differential  operator,  with  coefficients  independent  of  (R_,t),  while  Q, 
also  linear,  may  be  an  integro-differential  (or  global)  operator  and  is' 
dependent  on  position  CB,t).  (This  character  of  Q  arises  from  scattering 
element  doppler,  cf.  [36].)  For  the  moment  we  regard  Q,  and  hence  (2.8), 
as  deterministic. 

An  important  example  for  acoustic  propagation  in  underwater  media 
is  the  case  of  inhomogeneous,  absorptive  media  [2],  for  which  (2.8)  is 
specifically  [33],  [24] 


where 


(3.1) 


;<0)  e  <i  +  v 


X  <d_.  n  -  c  a2 
c23t2’  c2at2 


-  y  —  V*\ 
ox  x  9t 


(3.1a) 


with  t  =  t  • ( 1  +  y  [R,t])  generally.  When  x  =  0,  (3.1)  reduces  to 

A  UA  A  “•  U A 

the  more  familiar  "extended"  Helmholtz  equation 


wi  th 


■y  -  x  ci  ♦  s(R.tn  --j-jc.  = 
co  ~  at2 


’r 


(0) 


v2  1L 


(3-2), 


(3.2a) 
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Here 


(3.3) 


R  =  tx  +  ty  +  tz 
—x  ~y 

? 

and  rectangular  coordinates  are  assumed  throughout,  so  that  V*_V  =  V  , 
etc.  Specifically,  cQ  =  (constant)  wave-front  speed,  e  embodies  the 
effects  of  velocity  gradients,  internal  wave  phenomena,  and/or  (weak) 
local  turbulence;  while  xx  represents  the  effects  of  relaxation  absorption 
(due  to  I^SO^  and  other  salts). 

With  inhomogeneous  (but  still  deterministic)  media  conventional  tech¬ 
niques  of  solution  fail,  largely  because  the  medium  is  not  reciprocal, 
and  there  is  no  general  way  of  applying  and  evaluating  the  initial  con¬ 
ditions  over  the  various  volume  integrals  which  appear  in  the  development 
of  the  generalized  Huygen's  principle  (GHP)  when  Q  /  0  ([2]).  Moreover, 
since  these  media  are  space-time  variable  in  their  inhomogeneities,  e.g., 

A  /V 

Q  =  Q(R,t|  ...),  either  locally,  or  because  of  local  doppler  [36],  such 
media  do  not  support  space-  and  time-harmonic  solutions  [2],  nor  are 
the  standard  perturbational  and  variational  techniques  of  the  "classical" 

A  a 

approaches  [35]  general  ly  appl  icable,  particularly  when  Q  =  Q(_R,t|  ...) 
and  is  a  random  operator,  as  will  be  the  case  here.  Finally,  in  addition 
to  all  these  difficulties  is  the  usually  insurperable  one  of  bounded  media 
when  the  boundaries  themselves  are  random  and  moving. 


3.1  Random  Media 

Accordingly,  entirely  new  approaches  are  required,  for  manageable 
solutions  of  (2.8)  generally  (and  (3.1)  in  particular),  which  embody  the 
physical  realities  of  the  application  in  question  ([2],  [21]),  and  which, 
in  particular,  include  both  the  random  character  of  Q,  i.e.,  e(R,t)  in 
(3.1)  etc.,  and  the  random  doppler  effects  [36],  as  well .  These  involve 
incorporating  the  effects  of  boundaries  into  the  now  random  inhomogeneity 
operator  Q,  as  well  as  different  types  of  scatter  mechanisms,  e.g., 

0  =  %  *  «vol  *  ^Bottom’  etc"  with  5V0,  ■  QInh  +  %isc'  “here  ^Inh 
represents  the  continuously  (random)  variations  in  medium  properties, 

while  Qoisc)  embodies  the  effects  of  localized,  or  discrete  particles, 

bubbles,  etc.  A  typical  geometry  is  sketched  in  Fig.  3.1,  showing  the 

usual  variety  of  interactions. 
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Figure  3.1:  Schematic  scatter  channel  geometry  of  source  (T),  medium  (V), 
and  boundaries  (S,B),  with  receiver  (R),  D  =  unscattered,  or 
"direct"  path,  and  n  is  the  order  of  the  interactions. 
Sufficiently  high  frequencies  are  assumed  to  permit  illus¬ 
trative  ray  paths. 

-tor1 

For  these  linear  media,  with  L'  '  =  M  (which  also  includes  boundary 

and  initial  conditions),  it  is  easily  seen  that,  in  operator  form,  (2.8) 
becomes,  for  each  member  of  the  ensemble  of  propagation  equations: 

(3.4) 

where  la  =  a,  etc.  We  note  that  (3.4),  is  the  unscattered  field,  or 
field  of  the  purely  deterministic,  i.e.,  homogeneous  portion  of  the  medium 
here.  Moreover,  because  any  boundaries  are  included  in  Q,  M  is  replaced 
by  M  ,  so  that  now 

oo 

*H  =  ,  (3.4a) 

where  the  infinite  domain  operator  for  the  homogeneous  component  of  the 
medium  becomes 


'1  -  /,  dR'<i(R.t|R;.fU  )R.  .  (3  6a) 

—  GO  (  00  ) 


-T!(Y  )  =  }. 

*— s  v  o ,<»  k.  s  ^o,°° 


(3,5b) 


Y  =  Y  (s | p)  /  dt'e  st  (  ) .  , 

0,®  o,°°v  1  J  't 


(3.5c) 


k  =  1  k  +  t  k  +  t  k  -  2ttv 
— .  —xx  —y  y  —  2  z  — 


(3.5d) 


is  a  vector  wave  number,  with  a  vector  spatial  frequency,  cf.  (2.7). 

Here  g^  is  the  green's  function  of  the  corresponding  infinite,  homogeneous 
medium  (Q  =  0).  For  the  ocean  medium  supporting  propagation  according  to 
(3.1a),  we  have  specifically  the  (homogeneous)  operator  (kernels) 


'PF“(1+Tnxs)"  '"st‘ 
Lo 


4TTp(1+Toxs) 


(3.6a) 


tj  =  [k‘(l+t  s)  +  s2/c2]-1  ;  ps|R-R'U  kd  =  k-k.  (3.6b) 

^O,00  OX  '  '  oJ  '  —  — 

For  the  extended  Helmholtz  equation  (3.2),  (3.2a),  Eqs.  (3.6)  for  the 
homogeneous  part  reduce  to  the  simpler  relations 


-ps/c  -st' 


0,°°  4iTp 


■  (k2  -  4)*1; 

C^ 

0 


M  ( R ,  1 1  R '  ,  t ' )  = 


6 (t-t ' -p/ c  ) 


(3.7) 
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3.2  Field  Series  and  Feedback  Solutions 

The  ensemble  solution  of  (3.4)  is  now  given  formally  by  [2],  [21]): 

(3.8) 

where1!  is  now  the  random  field  renormalization  operator  (FRO),  when 
considered  over  all  member  equations  of  the  Langvin  equation.  As  long 

A 

as  ||T)  ||  <  oo  (for  each  member  equation),  (3.8)  supports  a  series  expansion: 

II  H  II  <  i;  (3-9) 

(or  a  =  -  l  (T]"1)(n+1)aH,  1  <  ||T]  ||<»).  (3.9a) 

n=0 

A  A 

(Usually  it  can  be  shown  that  ||T)  ||<1,  [2];  we  remark,  also,  that 

A 

and  Q  do  not  usually  commute,  but  commutation  can  always  be  achieved 
in  the  proper  transform  space,  usually  (k>,s)-space  [2].)  The  expansion 
(3.9),  or  (3.9a),  is  called  a  Perturbation  Theoretical  Series  Solution 
(PTSS),  with  the  convergence  condition  as  indicated  formally  above. 

Equations  (3.8),  (3.9)  are  exact  (ensemble)  solutions.  Eq.  (3.8), 
moreover,  is  called  the  Feedback  Operational  Solution  (FOS),  since  M 
may  be  interpreted  as  a  "feedforward"  and  Q,  a  "feedback"  operator, 
cf.  [2],  [11],  in  the  manner  of  Fig.  3.2.  In  fact,  the  PTSS  in  (3.9), 
expressed  as  a  series  of  loop  iterations  formally,  indicates  computa¬ 
tionally  how  exact  numerical  solutions  may  be  approached.  The  evaluation 
of  the  PTSS  may  then  be  viewed  as  both  a  problem  in  modern  control  and 
computational  theory,  extended  now  to  four  dimensions  and  governed  by 
partial  (rather  than  ordinary)  differential  equations  and  influenced 
by  a  radiation  condition  (on  M^) ,  as  well  as  boundary  conditions 
appl ied  locally  to  Q. 

However,  the  formal  solutions  (3.8),  (3.9)  are  not  very  useful  in 
their  present  form  since  f)  is  still  a  random  operator.  To  achieve 
useful  results  the  langevin  equations  must  be  converted  into  various 
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o(R,t) 


Figure  3.2  Feedback  Operational  Representation  (FOR)  for  the  propaga¬ 
tion  Eq.  (.3.8),  for  linear  inhomogeneous  media. 

deterministic  forms,  which  represent  the  various  moments  of  the  scattered 
field,  as  we  have  noted  above  in  Sec.  2.1.  This  is  done  in  the  manner  de¬ 
scribed  in  Appendix  A.l  and  leads  to  suitable  replacement  of  the  random 
Q  by  an  appropriate  deterministic  operator  »  etc.  for  the 

first-,  second-,  and  nigher-order  moments  of  the  field  x,  c f.  Secs.  A. 1-1. 
Thus,  we  get  a  Dyson  equation  for  <a>,  a  form  of  the  Bethe-Sal peter 
equation  for  <u  ot£>,  etc.  For  example,  the  FOR  of  Fig.  3.2  applies  now, 
with  a  replaced  by  <a>  and  Q  by  .  Then,  by  approximati  ng  o|d\  etc., 
the  various  terms  of  the  PTSS  can  be  evaluated,  although  this  is  not  an 
easy  matter. 

Finally,  we  observe  that  the  ensemble  solutions  (3.8),  (3.9)  may 
also  be  insightfully  described  by  equivalent  diagram  representations, 
which  are  also  particularly  useful  in  suggesting  appropriate,  simplifying 
approximations.  A  concise  summary  of  this  approach  is  provided  in  Appendix 
A.l,  some  results  of  which  we  use  below  in  the  text,  cf.  Sec.  5  ff. 

Related  methods,  involving  path  integrals  and  various  models  of  distributed 
ocean  inhomogeneities,  with  particular  attention  to  long  distance  sound 
transmission  in  the  ocean,  are  treated  in  the  recent  important  studies 
of  Flatte  et  al.  [46].  These  offer  a  variety  of  techniques  for  realizing 
the  inhomogeneity  operator  Q  and  its  associated  statistics  (see,  for 
example.  Chapter  10  of  [40].) 
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4.  Channel  Characterization  and  .Stati sties 

The  (random)  channel  operator (=  (i-Tj)'^M  ),  cf.  (3.8),  is 
next  conveniently  written 


t(N)  - 

i  +  i)  M  =  M  +  I  HT?{ 

i  •  24  s 

i-y\i 

I-TJ 

where  we  separate  the  homogeneous  (M  =  M  here)  from  the  scatter  operator, 
I,  cf.  (3.4)  et  seq.  We  note  that  for  ambient  fields,  in  addition  to  the 
desired  signal  source  (S-jn),  we  replace  the  source  function  -Gy,  in  the 
Langevin  equations  (2.8),  (3.4),  by  -Gy-G^,  where  now  is  a  localized, 
or  a  distributed,  ambient  "noise"  source  with  an  associated  field,  = 

M°o ( - G « ) .  In  any  case,  it  is  at  once  evident  that  the  medium  operators 
M,I  (•**7T| )  are  unchanged. 


4.1  Moments  of  the  Received  Field 

Mostly  for  reasons  of  analytical  complexity,  channel  characterization 
is  usually  formulated  in  terms  of  the  lower  moments,  e.g.,  the  mean  and  co- 
variance  operators  associated  with  I.  The  pertinent  moment  operators 
here  are,  directly  from  (4.1), 


=  M  +  <I>  ,  Kj  =  <1^1 2>  -  <ij>  <If>  ; 

<*lf)2>  =  m1m2  +  k/lf  <  +  <1^  , 


(4.2) 


and  with  higher  moments  similarly  determined.  These  moment  operators 
are,  of  course,  deterministic,  as  a  result  of  averaging.  (Mote  that 
j  :  I(Rj  2,t:i  etc-’  (R^T^),  (JR^)  bein9  (generally)  dif¬ 

ferent  points  in  the  field.) 

Expressions  for  the  correspondi ng  moments  of  the  received  wave  X(t) 
are  found  from 


(4.3) 
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Thus,  the  mean  received  waveform  and  various  second  moments  are  respectively 


<X>  =  R<a>  =  R(Moo  +  <!>){-Gt-Ga) 


(4.4) 


<X:X2>  =  RxR2  <0^  =  RjRg^GjGg,  G  =  GT  +  Gft,  (4.5) 


wi  th 


Kx(t1tg)  =  <X^2>  -  <Xj>  <X2>  =  R1R2(<a1a2>  -  <a.>  <^2>) 

=  RjR^itGjGg), 


(4.6) 


where  from  (3.8)  and  (3.9),  and  (4.1),  we  see  that  all  orders  of  moments 
of  the  inhomogeneity  operator  Q  appear  in  the  above. 

Equivalent  diagram  representations  for  the  scatter  operator  I  and 
field  covariance  operator  Kj,  (4.2),  are  given  in  Appendix  A.l. 


4.2  Remarks  on  Arrays 

The  array  operators  (2.4b)  and  (2.6),  which  appear  in  (4.3)  et  seq. 
for  the  received  wave,  can  be  rewritten  alternatively  as 

00 

R  5  /  dt  J  (  )h(n,t-T)pdn,  dV„  -  dn  .  (4.7) 

-0°  vR 


Here  VR  once  more  is  the  region  occupied  by  the  physical  array  itself  (in 
the  field  x)  and  n  locates  the  array  element  dn  with  respect  to  0D,  the 
coordinate  origin  of  (here)  the  receiving  array,  cf.  Fig.  4.1.  The  array 
is  a  space-time  linear  filter. 

In  practice,  most  physical  arrays  consist  of  an  assemblage  of  m  = 

1,  2,  ...,  M  discrete,  but  distributed  elements,  each  of  which  samples 
the  radiation  field,  a.  Thus,  the  array  operator  (4.7)  may  be  expressed 
in  more  detail  formally  as  a  vector  operator 


i*  !J'  L  (  >h,«(c-t-t>Rfc  ■  2*  'V  •  <4-3> 

m 
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which  defines  each  "component,"  R  .  The  relation  (4.8)  may  be  simplified 
further  to  the  idealized  but  still  very  useful  situation  of  "point" 
sensors,  e.g.,  maximum  sensor  dimension  (each  m)  is  always 

much  less  than  the  acoustic  wavelength.  Thus,  (4.8)  becomes 


Soint"  '  (6^m>  1 


dx  h  ( t-T  )  D} 
nr  'R 


h  05  (n-n  )) 
mn-R  L«*m' 


(4.9) 


all  m;  where 


Vn-R  =/  <  VT ’r*  (4.9a) 

-oo 

for  temporal  sampling  t  =  nAt;  n  =  1,  ...,  N,  as  well.  It  is  convenient 
to  use  j  =  m,n  =1,  . . . ,  J  =  MN  as  a  single  index  for  space-time  sampling, 
so  that  now 


i5  (iW  ■  <V  (4-I0) 

is  a  J-component  vector.  Consequently,  the  sampled,  received  field  X 
in  (4.3)  becomes  explicitly 


}  =  Ra  =  {/  btvVT)^^’!)^} 


+  n  . 


(4.11) 


Thus,  X  =  (X.  1  is  the  received  data  vector,  after  linear  space-time 

filtering  and  sampling.  Equation  (4.11)  is  the  basic  relation,  given 
the  field  <,  which  constitutes  the  received,  sampled  waveform  .  This, 
in  turn,  is  the  fundamental  input  to  the  "signal -processing"  systems 
noted  in  Fig.  2.1,  which  are  the  principal  topic  of  Part  II  following. 

For  narrow-band  signals  in  the  medium  we  can  provide  a  more  detailed 
structure  to  hm(  )D  above.  In  fact,  by  combining,  i.e.,  summing,  the 
m  (=1,  ...,  M)  outputs  of  each  receiving  sensor  with  appropriate  delays 
At  ,  we  can  form  a  beam  (in  the  far  field  of  the  incoming  signal),  e.g., 


X 


=  (lx  } 

m  J 


=  { y 

m 


v  -v  n ;  f  ) 
-o  ~oR  o' 


(  (t  -Ir 
m  n 


,)) 


(4.12) 
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Source 


) 


Figure  4.1 


Geometry  of  array  (in  VR),  and  m-^-- sensor,  with  beam  in 
direction  -i 

~OK 


where 

A<”>  = 

Arm  "  Tm  Tm  ilo  JoR^  Jrr/co  ’ 


in  which  v  -  =  t  nf  /c  ,  v  =  i  f  /c.  are  wave  numbers,  respectively 
~oR  ^-oR  o  o  *-o  ~'0  o  o 

associated  with  the  beam  direction  >tQR  ( =  ^oR/  IJSqr !  )  ar>d  the  direction 
of  the  field's  wave-front ,  vTQ  =  ^0/|£0l*  The  geometry  of  the  nJ—  - 
sensor,  the  beam,  and  signal  wavefront  are  shown  in  Fig.  4.1,  where 
specifically  the  unit  vectors  Jo’  JoR  are  given  by 
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In  (4.12),  |A^|  and  <}>Rm)  are  the  amplitude  and  phase  of  the  m^--  sensor's 
response,  at  the  center  frequency  (f  )  of  the  received  signal,  here  the 
field  a.  For  a  vertical  array  of  2M  +1 =  M  sensors,  for  example,  we  have 


r  =  l  mAS., 

'~-z 


( >  •  ♦  • »  m,  •  •  •  *  Mq j 


At  =  —  (sine  -  sine  ~), 
m  c  o  oR 


(4.14a) 


in  which  A £  is  the  sensor  spacing.  For  a  similar,  but  horizontal  array. 


r  =  i  mA£  and  .-.At  =  -  [cos<+>  cose  -  cos<i>  Dcose~] .  (4.14b) 

~m  '~x  m  c  L  vo  o  oR  oRJ 

o 

For  a  narrow-band  signal  field  t  is  usually  known,  or  findable. 

Thus,  the  "beam"  can  be  steered,  that  is,  pointed  in  the  direction  of  the 

desired  signal,  by  choosing _TqR  =  J  ,  or  setting  the  path  delays  xm  in 

Fig.  4.1  such  that  Axm  =  0,  cf.  (4.12a).  However,  for  an  ambient  or 

scattered  noise  field,  each  of  which  is  the  superposition  of  a  number  of 

possible  (undesired)  signal  fields  randomly  phased  with  respect  to  one 

another,  T  is  randomly  distributed  in  space.  There  is  consequently  no 

distinct  noise  field  wavefront  or  wavefront  direction,  and  f  must  be 

~'0 

averaged  over  in  suitable  fashion.  [See  Secs.  8.2  and  Sec.  9  of  [4],  for 
example;  also,  (6.10b)  and  Sec.  7.3  ff.].  This  randomness  of  l  in  the 
case  of  noise  fields  usually  leads  to  spatial  non-uniformity  of  the  inter¬ 
fering  field,  which  can  often  be  exploited  to  increase  the  processing 
gain  in  signal  reception,  by  effectively  increasing  the  number  of  inde¬ 
pendent  spatial  noise  samples  vis-§-vis  the  (nearly)  uniform  desired 
signal  field,  cf.  Sec.  10  of  [4]. 


5_.  0cean^Vol ume  and  Surface  Channels: 

The  Received  Field  X 

Both  the  important  cases  of  weak  volume  scatter  in  the  ocean  and 
ocean  surface  scatter  are  well  approximated  by  at  most  the  scattering 
element  interacting  with  the  homogeneous  field,  the  former  because 
multiple  scatter  effects  are  usually  quite  ignorable,  the  latter  from 
the  physical  geometry  of  the  wave  surface,  which  likewise  discourages 
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multiple  scatter  except  within  a  typical  wave  crest- to-crest  domain,  cf. 
Fig.  3.1,  (n=l).  (Similar  considerations  usually  apply  for  ocean 
bottoms,  as  well,  cf.  Fig.  3.1.  We  shall,  however,  consider  only  the 
modes  M=S,V,  here),  observing  also  that  o(30-40  db)  they  are  effec¬ 
tively  independent  ([21] ,  thesis):  there  is  negligible  coupling  in  the 
ocean  between  surface  (bottom)  and  volume  scatter.) 

Accordingly,  a  lst-Born  approximation  suffices  (as  indicated  by  the 
dotted  lines  in  (A. Id),  (A. 2a),  and  we  write  for  the  received  field,  X, 
from  (A. Id): 


The  mean  wave,  <X>,  is,  with  <Q^  =  •  ,  cf.  (A. 2a): 

<X>  =  R<o>  =  R  ( - >(f)+  ->#■■  '■>£)  )  =  j  RaH+RM<QM>aH  | 

+  aT  S  >  .  (5.2) 

(XH  +  XI  1 

The  second-moment  of  X,  (5.1)  in  (4.5),  is  found  to  be,  on  averaging 
over  the  phases  (4>)  of  the  input  signal. 


<<X1V>,»  ’  “  <S1fe2['1  *  V2 


•  *  j  <q'k)q'k)>] 


(5.3a) 


XHlaH2^>s 


<<X1X2>>?  =  <R1R2(| 


k  =  0  k  =  1  (k  >  2) 

>0  — ►  X>  +  — > 

11  22  11  22  11  22  "12  ^1 

:oherent  +  incoherent  multiple- 

scatter:  incoher. 


11  22 


"homog. "  \  "classical " :  coherent  +  incoherent 


unseat . 


coh.  +  incoherent  (0) 


(5.3b) 
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where 


O  =  q  =  Q  -  <0>R  =<— 0-0  =*{v(0)  =  l  «v^k^>R+Av^kb}; 

K  k=l 

Q(0)  =  <Q>r  =  o,  .*•  <Q(0)>S  =  <Q>  =  •  =  <0>R. 

space  av.xe  5  =  ^lV  etc-^^=  <  >  =  <  >r,s  ’ 

(5.4) 


via  (A. 2-6).  Note  that  the  (total)  average,  <  >,  as  in  Appendix 
A.l  et  seq.  ,  involves  both  an  average  over  radiation  events  and  an 
average  over  random  spatial  (S)  and  parameter  values  (e ) ,  embodied  in 
<><.,  cf.  (5.4)  and,  particularly,  Appendix  A. 2. 

^  A 

It  is  the  Decomposition  Principle,  embodied  in  Q,,,  (M  =  S,V,B),  Eq. 

(A. 2-5)  et  seq.,  which  allows  us  to  resolve  the  inhomogeneity 
operator  into  sets  of  distinct  and  (statistically)  independent  entities 
as  specifically  exhibited  in  Eqs.  (5.3).  Moreover,  and  most  important, 
the  entities  (k  =  0,  k=l,  k>2)  in  (5.3)  above  (and  in  (5. 1) ,  with  the 
help  of  (5.4))  have  an  explicit  physical  interpretation:  Thus,  the 
terms  k  =  0  contain  both  the  coherent  and  incoherent  radiation  contribu¬ 
tions,  from  all  orders  {k>  1)  of  radiation  interactions,  e.g.,  single- 

(kl 

and  multiple  scattering  (k^2),  for  example.  Consequently,  if  vv  ' 
represents  the  density  of  (illuminated)  k-coupled  scattering  elements 
(k>l),  then  <vv  '>R  is  the  coherent  radiative  contribution,  where 
<  >R  is  the  average  over  the  radiation  events  associated  with  the  ensemble 
of  potential  (re-)radiating  sources.  In  a  similar  way 


Av 


(k) 


(5.5a) 


is  the  fluctuation  in  the  density  of  k-coupled  radiating  elements,  and 
is  always  associated  with  incoherent  radiation.  Accordingly,  we  have 


v 


(0) 


l  (<v(k)>D+Av(k) 
k=l  R 


)  = 


oo 


I  V 

k=l 


(k) 


(5.5b) 
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showing  the  basically  coherent  and  incoherent  elements  of  the  scattering 
or  source  region  containing  these  various  k>  1  types  of  elements.  The 

(average)  density  of  coherent  radiators  is  <v^>D  =  y<v^)>  since 

(k)  K  K 

<Av '  ;>R  =  0.  The  number  of  "radiation  events"  occurring  in  a  small 

region  dA  of  the  illuminated  or  emitting  domain  A  is  dN,  so  that 

v  =  dN/dA  =  v{0).  A  taxonomy  and  interpretation  for  the  dN,  similar  to 

(5.5a,b),  can  be  constructed:  One  has  now 

dN^  5  l  (<dN^kVD+dn^k^)  =  l  dt^k^.  (5.5c) 

k=l  K  k=l  . 


The  relations  (5.6)  et  seq.  obviously  apply  if  one  replaces  v  etc.  by 
dN;  see  also  Sec.  6.1  following. 

Since  the  medium  in  question  is  linear,  the  scattering  or  inhomo- 

-  f  k) 

geneity  operator,  Q,  is  a  linear  functional  of  the  i.e.,  of 

v(0>,  e.g.,  Q  =<£(v^),  thus  containing  all  types  of  radiation  inter¬ 
action  effects.  Moments  of  v^,  or  more  generally,  of  (L,  which  appear 
in  (5.2),  (5.3)  above,  are  readily  indicated:  we  have 

<Q>  =  £{<v^>}  ,  where  <\Z°^>  =  «v^°^>R>s.  (5.6) 

This  is  different  from  zero,  in  scattering,  if  specular  reflection 
geometry  (both  on  a  surface  and  in  a  volume)  is  available.  The  second- 
order  moments  are  instructive: 

<Q1Q2>  =2«v[°40)>}  ,  (5.7a) 


where  now 


<v 


(0)  (0) 

1  v 


>  =  «Av 


(0) 

1 


*v2° Vs 


«v<0)>R<v2 


(0) 


(5.7b) 


Equation  (5.7b)  contains  both  incoherent  (~Av^)  and  coherent  terms 
(~<Av^>R),  of  all  orders  (k^l)  of  coupled  scattering  or  radiating 
elements.  The  associated  covariance  of  Q  is 
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(5.7c) 


^«vj0)v^0)>  -  <v[0)>  <v^0)>} 


which  embodies  wholly  incoherent  radiation. 

We  now  distinguish  two  principal  classes  of  radiation  models  from 

which  scattered  and  ambient  acoustic  (and  electromagnetic)  noise  fields 

can  be  constructed.  These  are  (1):  the  "classical"  cases,  where  the  number 

of  "radiation  events"  (e.g.,  reradiations  or  ambient  emissions)  are  the 

(k)  (k) 

same  for  each  member  of  the  ensemble,  so  that  <vv  '>D  =  ,  and 

fkl  K  (k) 

•  ■  Av'  '  =  0,  all  k>0.  And  (2):  the  quantized  cases,  where  <vv  \>D, 

(kl  K 

and  •'•A v  ,  can  be  nonvanishing,  reflecting  the  fact  that  the  number  of 

radiation  events  is  variable  over  the  ensemble.  The  latter  is  typical 

of  the  so-called  poisson  radiation  models  [11],  developed  recently  for 

both  scattering  [28]-[32],  [40]  and  ambient  noise  [3]-[6],  [52]. 

For  the  poisson  scattering  cases  there  is  the  so-called  FOM  (Faure 

[28],  01‘shevskii  [29],  Middleton  [31])  model,  based  on  an  independent 

poisson  point-process,  where 

(FOM):  v(0>  -  v(l>  ;  =  0;  k  >  2,  (5.8) 


namely,  all  multiple  scatter  is  neglected.  For  this  model  the  coherent 
and  incoherent  contributions  are  found  to  be  linear  functionals  of 


<v^>R  =  <v^>R  ,  or  <v^>  =  <v^>  i  ; 


(5.9a) 


(FOM):  K  ^  =  <v(°)v(°)>  -  <v[0)>  <v^> 


<o(1)^2-il»S  *  -  <V1U>  <u21)>' 

(5.9b) 


where  z. ,  z0  embody  spatial,  temporal  and  parametric  coordinates  and 
(1) 

ov  y(z)  is  the  average  density  of  scattering  in  the  illuminated  region  A. 
The  FOM  model  is  quasi -phenomenological  and  avoids  dealing  explicitly 
with  the  boundary  conditions,  which  are  now  contained  in  a  (linear)  ad  hoc 
response  function  Y^t^),  with  "cross-section"  yq,  developed  and 
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empirically  adjusted  for  the  problem  at  hand  [43]-[45]. 

On  the  other  hand,  in  the  classical  approach,  the  density  fluctuations 
of  radiation  events  are  always  zero,  as  noted  above,  so  that  (5.7b), 
<vi°M0)>  =  is  an  identity,  and  (5.7c)  describes  the  covariance 

( Qr  /s  l  ^ 

of  vv  '  ,  or  Q.  Again,  in  the  usual  applications,  multiple  scatter  is 
(kl 

neglected  (vv  1  =  0,  k>2).  One  has,  here  for  surface  scatter  cases, 
specifically 


Vvl0)>R  * 


V 


(0)  i  Y  vd)  (?  -?  ) 

'o  c  ~g  ~T  •^R 


(5.10) 


under  a  Kirchoff  approximation,  [27],  [28],  where  R,S  are  respectively 

the  (average)  plane-wave  reflection  coefficient  and  shadowing  function, 

oj  (=2Trf)  is  the  angular  frequency  of  the  incident  radiation,  cQ  =  phase 

speed  of  sound  in  water,  nn  =  (i +i,  4. -i ,)//l+C+C  is  the  normal  to 

x  ~y  y  x  y 

the  wave  surface  ^(jr,t),  with  z,  =  3 ^/9x,  etc.,  and  j,,  are  unit 
vectors  from  source  and  receiver  to  the  scattering  element  on  z,.  See 
[25]- [27 3 ,  [37]-[39]  as  well  as  [47]-[^l],  for  detailed  development  of 
the  classical  theory.  See  Sec.  6  following  and  Appendix  A. 2  for 
further  discussion  of  these  poisson  field  models. 


6_^^Poisson  Statistical  Models:  Non-Gauss ian 
Noi_se_Fjelds  ajid  Receij/ed  Processes,  X 

Because  of  the  effectively  discrete  nature  of  the  inhomogeneity 
models  and  the  analogous  structure  of  the  ambient  field  sources,  the 
fundamental  statistics  of  both  the  fields  u(R_,t)  generated  in  the  medium 
and  the  received  waves,  X(t),  are  poissonian.  This  is  directly  the  case 
for  the  ambient  noise  and  weakly  scattered  field  components,  which  consti¬ 
tute  much  of  the  interference  in  practical  operation.  In  many  cases, 
i.e.,  where  there  may  be  a  few  relatively  strong  sources  or  scatterers 
outstanding  from  the  majority,  the  resultant  field  and  received  process 
X  are  highly  nongaussian,  while  the  weaker  background  is,  as  expected, 
well  described  by  a  limiting  gaussian  process.  We  note  some  of  the  im¬ 
portant  results  below. 

6.1  Poisson  Helds  and  Statistics 

A  (scalar)  poisson  field  <(^,t),  here  an  underwater  acoustic  field 
generated  by  potentially  many  independent  sources,  may  be  represented  by 
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a(R,t) 


f  *(R,t|z)dN(Z(z)), 


(6.1) 


where  ,  in  which  ^zR  represent  radiation  epochs,  ^  are 

spatial  variables,  and  0^  are  structural  or  waveform  parameters; 

Z  =JRx<Zsx<Ze  is  the  associated  composite  space  of  these  quantities;^ 
is  the  wave  field  for  a  single  emitting  source.  Here  dN  is  a  counting 
functional,  and  in  general,  dfl,  z^,  0  are  random  (process  and) 
variables.  This  counting  functional  is  represented  by 


dN(Z) 


I  6(x-x  )6(e-e  )6(t-t  ) , 

-  -w*  w^n  w  >~n  n 


p=i 


(6.2) 


for  P  discrete  sources  potentially  emitting,  or  reradiati ng, at  times 
t  in  region  A,  where 


dN(Z)  =  dN(A=p;t  ;6)  =  (dZ)Pe"^/p!,  p>  0  (6.2a) 

~~  p  ~~  - - 


with  dZ  =  osdAde '  dt'w^ej  )w^(t' ) ,  etc.,  in  which  t'  is  a  random  emission 
time  or  epoch.  For  poisson  fields,  dN  is  a  poisson  process  obeying  (5.2a) 
such  that  for  dn  =  dN  -  <dft>. 


e.g. , 


Q-l 


<dn^. .  .dn^  =  p(z^1)dz1  n^6(Zj+1-z^)dZj+1;  <dn^  =0;  (6.3) 


<dn1dn2>  =  p  (z^i^-^dz^. 


(6.3 a) 


where  p(>  0)  is  the  "process  density,"  cf.  (A. 2-12).  Thus,  for  the  mean 
field  <o>  and  field  covariance  we  have,  on  applying  (6.3)  to  (6.1), 
the  familiar  relations 


<a(R,t)>  =  /  ^(R,t|z)pU)dz; 


K((Ei’t1’^2’t2^  =  /^Rrti!z.)^(5.2>t2ll)‘1(z)dz, 


(6.4a) 

(6.4b) 
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with  analogous  results  for  the  higher  moments.  Since  X  =  Ra,  cf.  (2.4b), 
(5.1),  the  corresponding  results  for  the  received  process  X  are 


<X(t)>  =  /  L(t|z)p(zjdz,  L  =  R& 

KX(tl’V  =  / 

The  Q^-- order  characteristic  function  (c.f.)  is  readily  found 
to  be  ([10]  ,  [31]  ,  Part  I) 

1^9  Tq(iC^i*--5lCQ|o!^»-**'>ct^  or  X^,...,Xq)  - 


J  (6.5a) 

]  (6.5b) 


Q 

<  / p(zj  ,z)[exp{  i  l  or  L£)}  -1]  dz>z , ,  (6.6) 

2,=  1 


where  z  =  (z',z)  and  z  =  t/Tc  is  a  normalized  time,  with  i£.  =  JC(R.  ,t,  I  z) , 

j  l  ““i  j.  “■ 

L1  =  1(1^^),  etc.  [The  limiting  gaussian  cases  follow  formally  from 

(6.6)  by  expanding  the  exponent  through  )  and  dropping  higher- 

th  K  ^ 

order  terms.]  The  associated  Q - order  probability  density  function 

(pdf)  is  formally,  with  y  =  <£or  L: 


wj(yr-  *  •  ’-V 


/  exp  v  >Fn(U|y) 


(2tt)C 


dC 


dV  ■  ,d^Q’ 

(6.7) 


Evaluating  (6.7)  based  on  (6.6)  is  generally  nontrivial.  Only  the  cases 
n  =  1,2  appear  analytically  manageable,  but  fortunately  w^(y)  is  suffi¬ 
cient  to  give  us  useful  results  in  signal  processing  appl ications ,  as 
we  shall  observe  in  Part  II  following.  [The  case  Q  =  2,  e.g.,  w^y^y.,) 
allows  us  to  extend  the  "classical"  theory  of  signals  and  noise  through 
(zero-memory)  nonlinear  devices  to  include  specific  nongaussian  noise 
processes  and  fields,  as  recent  work  has  shown  [ 10] . ]  We  remark  that 
our  ability  to  construct  the  field  covariance  (6.4b)  allows  us  to  estimate 
the  degree  of  spatial  processing  gain  (1  ^  M1  s?  M)  attainable  when  the 


etc . 
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field  is  nonuniform  in  space.  Here  M'  is  the  estimated  number  (not  neces¬ 
sarily  integral)  of  effectively  independent  sensors,  out  of  a  total  M,  cf. 
comments  in  Sec.  1  above,  and  Refs.  [3]-  [5];  also  Sec.  8.3  ff. 

As  a  specific  example,  for  the  many  "far-field"  situations  where  it 
is  meaningful  to  talk  of  beams  and  beam  patterns,  cf.  Sec.  4.2,  the  basic 
(complex)  field  X  can  usually  be  described  (in  the  far-zone)  by  an  expres¬ 
sion  of  the  form 


/  s  (f(z  )aTU|f)e 

-CO 


iu)Pd(t-e-|R0-R|/c0) 


df, 


(6.8) 


where  now  X  in  (6.1),  (6.4),  (6.6)  is  X  =  Re  X.  Here  the  typical  source  is 

treated  as  essentially  a  point  source,  with  a  non-uniform  beam  pattern  Q.y, 

Fig.  4.1,  and  X  =  |R  -R|/c„,  with  X  =  |R  -r  j/c  the  distance  of  the  m— - 
3  l'*-o  — 1  o  m  'v-o  x»m  o 

array  element  from  the  typical  source,  and  y>0.  Here,  also. 


„  /  00  00  i (t-x)  ,  -icot 

ST(f|z)  =  [f  /  a(t-x)sj(x)e  dxje  dt,  u  =  2*f,  (6.9) 

-CO  -CO 

represents  the  source  waveform  s-j-  with  a  fading  mechanism,  ae  ;  ud 
(s  l+ed)  =  doppler  coefficient  in  (6.8),  while  z  (=T-t)  is  a  (normalized) 
time,  cf.  (6.6). 

From  the  above  we  can  show,  for  instance,  that  the  covariance  of  the 
interference  output  of  the  receiving  array,  i.e.,  X(t)  =  Ra,  (4.12),  be¬ 
comes,  for  these  sources  narrow-band  about  frequency  f 


H(t)R  2  B  A(A  or  B)<02y>Re'-Ma^T^Dl^o;  *°(x)  <  ^R^o'^oR’VI  >i 


~o 


(6.10) 


Vi  1UJ  pT 

where*M  (t)  =  complex  covariance  of  the  fading,  ae  ;  D.(lj  r)  =  <e  0  > 

i  0 


is  the  doppler  "smear"  factor,  and 


£d 


~o 


y  W  /Ao ^o'-loR^  * ^ 

L  <e  >  , 


m  ,m 


o'  ro 


(6.10a) 


*  A 

f,0  is  the  typical 


interfering  source  covariance  (of  s() 
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represents  the  averaged  beam  pattern,  over  all  the  different  angles  (© 0 »4> 0') 

of  arrival  of  the  various  interfering  source  wavefronts;  Am  =  (complex) 

element  weighting,  A  =  normalization  factor,  for  example,  for  an  equally 

2 

weighted  vertical  array  (4.14a)  when  wi(60’4>0)  =  1/tt  ,  -t/2 <  0Q^ it/2; 
0<4>o<tt,  i.e.,  (e0 )  are  each  uniformly  distributed,  for  the  mm'  element 
in  (6.10a)  we  get 


-  e2TriA^,sine°R 

K  K  1 

w*"0 


Jo(2tZW>'  4r«n'  =  o' 


(6.10b) 


where  A 2.  is  the  interelement  spacing  and  \Q  (=  2ttcq/wo).  Generally,  the 
random  angles  of  arrival  "smear"  the  beams.  Moreover,  Eq.  (6.10b)  in 
(6.10a),  (6.10)  shows  at  once  the  (statistical)  nonuniformity  of  these 
interference  fields,  as  seen  over  a  typical  array,  when  the  array  is 
large  enough,  e.g.,  2MA2./Xq  =  l/Xq  =  0(1  or  more)  here,  cf.  J0(27TAmm./^0) 
in  (6.10c),  etc.:  the  field  covariance  as  sampled  by  the  array  elements  is 
clearly  quite  nonuniform.  [We  shall  return  to  this  point  in  Sec.  7.3  ff.] 

6.2  Non-Gaussian  Process  PDF's 

As  we  have  emphasized  in  previous  work  [4],  [6]-[8],  [10],  [13], 

[15],  [16],  [52],  involving  EMI  interference,  ana  analogously  here  for 
underwater  acoustic  applications,  most  noise  environments  can  be  char¬ 
acterized  in  three  main  classes  of  nongaussian  process  or  field:*  Class  A, 
or  "coherent"  noise,  producing  negligible  transients  in  a  typical  receiver; 
Class  B,  "impulsive,"  i.e.,  broadband  vis-a-vis  the  receiver,  generating 
transients;  and  Class  C  =  Class  A  +  Class  B,  in  various  combinations. 
Typically,  Class  A--for  example,  reverberation,  ship  noise,  etc. --or 
Class  B  (biologic  (shrimp),  ambient  ice,  microseisms,  impulsive  ship  noise, 
etc.),  predominate  in  most  instances. 


*An  important  exception  is  the  case  of  one  (or  a  small  number)  of  inter¬ 
fering  signal  sources  of  known  emission  behavior,  whose  deleterious  effects 
can  often  be  removed  by  selective  temporal  (=  frequency)  or  spatial  deletion. 


29 


For  our  subsequent  detection  and  estimation  algorithms  (Part  II  ff.), 
we  shall  need  a  more  sophisticated  statistic  than  just  the  process  or 
field  covariance  functions,  although  these  are  always  useful.  Specifically, 
we  shall  need  the  (first-order)  probability  density  of  the  (total)  inter¬ 
ference,  X(t),  w^(x|Hq),  as  obtained  from  (6.7),  since  we  are  going  to 
sample  in  such  a  way  as  to  produce  independent  noise  samples,  cf.  Sec.  7 
ff.  For  this  we  set  Q=1  in  (6.6)  and  carry  out  the  resulting  evaluations, 

along  the  lines  described  in  [6],  [8].  The  result,  finally,  for  Class  A 

t  h 

interference  is  found  to  be,  typically,  at  the  m— -array  element  now 


'A+G 


-a 


(m)2r2/2-A, 


AAe 


(fTI  )  /  o 

-  ,  -*2A  U 


[i  +  f[m)(-;2)A], 


(6.11) 


where  =  mean  intensity  of  the  gaussian  component,  =  mean 

intensity  of  the  (usually  much  stronger)  nongaussian  interference,  both  at 
the  nr—  array  element.  Here  because  of  the  assumed  (local)  homogeneity 
of  the  noise  field,  =  Oq,  2^  =  A^  =  overlap  index  or  "usage” 

parameter  [7],  which  is  typically  o(O.l-l)  for  nongaussian  underwater 
acoustic  interference.  The  quantity  is  a  correcting  factor,  which 

can  noticeably  modify  the  "tails"  of  the  associated  pdf,  cf.  (6.11), 
(6.12),  when  the  potentially  interfering  sources  are  widely  distributed 
in  space  ("Quasi-canonical "  cases,  [8]). 

When  these  sources  are  not  so  widely  distributed,  =  0,  and  we  have 
the  more  familiar  "canonical"  and  approximately  canonical  Class  A  cases. 
For  these  latter  the  pdf  associated  with  (6.11)  is 


(6.12) 


where  r*  •  and  x,,,  „  • 

is  the  normalized  input  data  sample,  at  array  element  m  and  at  time  t  . 
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(For  the  "anatomy"  of  the  noise  parameters  r^,  i^,  see  [7],  [8].) 
Class  B  interference  is  similarly  developed:  we  find  that  [6] 


wl(xm,nlHoW 


r( 


pg+1 

2 


i  x 

c  /  ga.  1.  m,n, 
i  i  2  ’  2’  r.  ' 


(6.13) 


Herea=  is  a  parameter  based  on  source  distribution  (~\  u )  and  propaga¬ 
tion  law  (~a~y),  while  ft  is  a  normalization  factor  [chosen  to  cut  off  the 
tails  of  (6.13)  at  some  very  small  residual  probability  and  insure  that 
/Wjdx  =  1:  (6.13)  does  not  support  a  finite  intensity,  e.g.,  x^-*  unless 
the  pdf  is  suitably  truncated.] 

Finally,  we  emphasize  again  the  canonical  nature  of  these  nongaussian 
models:  the  forms  of  the  pdf's  are  independent  of  the  particular  noise 
mechanisms  involved.  Of  course,  for  specific  applications,  we  must  esti¬ 
mate  the  relevant  parameter  values  and  insert  them  appropriately  in  the 
chosen  algorithm.  Thus,  we  may  say  that  signal  processors  based  on  these 
noise  models  are  parametrical 1y  adaptive,  beased  on  the  above  statistical - 
physical  models  of  the  acoustic  noise  or  interference  environment. 


6 .3  Summary  Remarks  (Part  I) 

In  the  preceding  sections  (Part  I)  we  have  outlined  a  formal  apparatus, 
illustrated  by  some  simple  examples  and  specific  results.  This  apparatus 
provides  a  general  approach  for  obta i ning  the  needed  noise,  signal  fields, 
ana  received  process  models,  which  constitute  the  inputs  to  our  sub¬ 
sequent  signal  processors,  namely,  signal  detectors  and  estimators.  Key 
features  of  these  inputs  which  must  be  taken  into  account  are  (i)  the 
generally  nongaussian  character  of  the  accompanying  noise,  or  additional 
acoustic  interference  from  localized  sources;  (ii)  the  spatial,  as  well  as 
temporal  character  of  the  noise  and  signals;  (iii)  the  salient  features  of 
propagation  in  a  typically  inhomogeneous  medium  which  includes  discrete, 
continuous,  and  distributed  inhomogeneities;  and  (iv)  the  spatial  coupling 
to  the  medium  itself.  The  controlling  factors,  as  always,  are  1  1timately 
physical.  These,  in  turn,  may  be  expected  to  impact  strongly  upon  the 
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design  and  operation  of  the  subsequent  signal  processing,  as  we  shall  see 
in  Part  II  following.  Finally,  but  not  least,  we  stress  again  the  direct 
analogy  between  the  underwater  acoustic  field  and  process  models  outlined 
here  and  corresponding  EMI  models  developed  for  telecommunication  appli¬ 
cations:  the  formal  structures,  and  ambient  noise  models  (neglecting 
multiple  scatter)  are  the  same,  as  are  the  basic  processing  concepts, 
so  that  results  from  the  latter  discipline  can  be  translated  to  the 
former,  with  appropriate  attention  to  the  governing  physics.  In  any 
case,  our  approach  emphasizes  the  concept  of  the  generalized  channel, 
in  the  manner  of  Fig.  2.1  above,  where  in  addition  to  propagation  ques¬ 
tions,  we  include  the  necessarily  coupled  signal  processing  ones,  cen¬ 
tered  on  the  critical  limiting  threshold  cases. 

Part  II.  Acoustic  Threshold  Signal  Processing 
7.  Threshold  Signal  Detection 

V"  w  '  .1  -II  I  “"‘■y  .rw>  —  

We  are  now  ready  to  outline  the  aforementioned  canonical  approaches 
to  signal  processing,  when  the  desired  signal  is  weak  compared  to  the 
noise  background.  Here,  of  course,  the  desired  signal  (and  undesired 
noise)  are  the  received,  scattered  and/or  ambient  acoustic  fields,  now 
sampled  by  our  distributed  arrays  (cf.  Sec.  4.2).  Our  principal  tasks 
are  to  detect  a  desired  signal,  and  then,  having  determined  its  existence 
in  the  accompanying  noise,  to  extract  some  one  or  more  desired  attributes 
of  the  signal,  namely,  to  perform  so-called  "signal  estimation." 

As  before  [7],  [15],  [16],  we  focus  our  attention  primarily  on  the 
"on-off"  threshold,  or  weak-signal  detection  situation,  involving  now 
nonuniform  interference  fields  (I),  where  the  decision  process  is  the 
hypothesis  test:  H^.Sel  vs.  Hq : I  alone.  A  canonically  optimum  theory  is 
again  possible:  canonical  in  the  sense  that  the  formal  results  for  de¬ 
tector  structures  and  performance  are  independent  of  any  detailed  physical 
mechanism,  as  well  as  of  specific  signal  waveforms  and  noise  statistics 
[16].  Although  the  results  are  no  longer  optimum  for  stronger  signals, 
they  are  generally  absolutely  better  than  for  weak  signals,  which  in 
turn  is  normally  satisfactory  for  applications.  Moreover,  an  optimal 
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theory  suggests  reasonable  suboptimum  procedures  which  can  have  the  advan¬ 
tage  of  structural  and  operational  simplicity,  as  well  as  performance 
close  to  optimum  [15],  [16]. 

Three  modes  of  detection  are  noted:  (i)  coherent,  when  signal  epoch 
is  known  precisely  at  the  receiver;  (ii)  incoherent,  when  signal  epoch 
is  unknown;  and  (iii)  compos i te ,  in  which  a  sum  of  modes  (i)  and  (ii)  is 
employed,  to  take  advantage  of  any  signal  epoch  information  available  [16]. 
Since  weak  signals  are  postulated,  large  processing  gains  (here  space-time- 
bandwidth  products  J)  are  needed  to  achieve  acceptably  small  probabilities 
of  decision  error.  This  ensures  that  the  detection  algorithm  in  question 
will  be  asymptotically  optimum  (AO)  (as  or  J  >>  1  practically), 

and  normally  distributed  under  Hg,  provided  a  suitable  bias  term 
(independent  of  the  data  and  determined  a  priori)  is  employed.  (The  im¬ 
portance  of  the  proper  bias,  generally,  is  stressed  in  [16].)  These  re¬ 
sults  enable  us  to  calculate  performance  directly  in  canonical  form. 

For  specific  applications  we  must,  of  course,  "calibrate"  our  canon¬ 
ical  algorithms  to  the  particular  interference  environment.  This  requires 
establishing:  (1)  the  Class  (A,  B,  or  C)  of  interference,  and  (2)  esti¬ 
mating  the  relevant  parameters  of  the  associated  probability  densities 
Wj(x)^,  Wj(x)B,  etc.,  [15],  [54],  This  can  be  done  if  we  can  construct 
the  EMI  scenario,  for  example  [4],  [8],  or  directly  by  empirical  observa¬ 
tion  [6],  [7].  Thus,  these  (detection)  systems  are  adaptive  and  include 
adaptive  beam-forming  here.  Their  often  considerable  improvement  over  con¬ 
ventional  (i.e.,  "matched-filter"  or  correlation)  detectors  stems  from  this 
essential  feature,  [7],  [16]. 

In  this  necessarily  concise  overview  we  confine  our  effort  to  an  out¬ 
line  of  some  of  the  principal  results  of  the  extensions  of  our  earlier 
work  to  include  spatial  processing  [3],  [4].  The  present  section  is 
organized  as  follows:  Sections  7. 1-7.4  provide  examples  of  general  and 
specific  space-time  processing  algorithms,  the  latter  for  independent 
sampling;  [Section  6  earlier  presents  some  results  for  nongaussian  field 
models].  Section  8  ff.  treats  performance,  in  both  optimum  and  suboptimum 
cases.  Section  8.3  completes  our  detection  treatment  here,  with  a  brief 
discussion  of  the  principal  assumptions  and  conditions. 


33 


7.1  Canonical  Space-Time  Threshold  Algorithms 

To  obtain  the  desired  optimum  (binary)  decision  algorithms  we  form 
first,  as  usual  [20],  [54]»  the  generalized  likelihood  ratio,  now  extended 
to  include  spatial  as  well  as  temporal  sampling,  e.g., 


£*  e  log  AT(xje)  =  g^xje)  +  t* ( X )  e ) , 


(7.1) 


where  J  =  MN,  with  N  =  number  of  temporal  samples  of  the  received  data 
=  (x....x._  ...x,),  at  each  of  the  M  sensors,  e.g.,  x  =  (x  , 

xmn_j  xmN),  with  J  =  1».  • .  ,mn,. . .  ,J=MN,  alternatively,  over  all  sensors. 
Here  0  denotes  an  input  signal -to-noi se  ratio;  g*,  t*  are  respectively 
the  desired  algorithm  and  a  remainder  series,  which  vanishes  prob.  1 
under  Hq,  if  certain  conditions  are  satisfied  [cf.  Sec.  7.4  f f .] . 

For  these  threshold  cases  log  Aj  is  expanded  aboutQ  =  0,  with  one 
or  two  terms  in  the  data  (xj  retained,  depending  on  the  mode  of  reception. 
The  required  bias  term,  B*(q),  is  obtained  from  the  average  over  { x}  with 
respect  to  the  null  hypothesis  HQ  of  the  next  nonvanishing  term  in  the  gj- 
expansion.  Thus,  we  write  formally 


(7.2) 


where  for  the  coherent  and  incoherent  detection  modes  we  have 


9J^coh  "  0F1^J  +  Bj(9  )Coh;  BJ-coh  =  2T<F2^d>Hn  +  log  M; 


(7.3a) 


9J^inc  =  2!  F2^J  +  BJ^9  ^inc;  BJ-inc  =  ;^>HQ  +  log  U' 


(7.3b) 


The  specific  structure  of  F^,  F^  depends,  data-wise,  on  the  basic  pdf 
Wj(x|0=O)  of  the  interference  alone,  and  on  the  desired  signal  waveform 
structure.  Here  u  =  p/q ,  p+q  =  1,  (p,q>0),  where  p,q  are  respectively  the 
a  priori  probabilities  that  the  data  sample  x,  contains  (or  does  not  contain) 
a  signal.  The  decision  process  (5)  is 
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decide  if  g^>log  K  or  decide  Hq  if  g^<log  K  (7.3c) 

where  K  is  a  threshold  in  the  usual  way. 

Both  because  of  much  needed  technical  simplicity,  and  the  fact  that 
it  can  be  shown  [4]  that  comparatively  little  further  improvement  in  per¬ 
formance  is  obtained  by  fully  correlated  sampling  vis-a-vis  (non-sparse) 
independent  (space  and  time)  sampling,  we  focus  henceforth  on  this  latter 
case.  Moreover,  although  the  theory  has  been  developed  for  the  general 
case  of  nonstationary  and  inhomogeneous  fields  [4],  we  shall  describe 
results  here  only  for  the  (often)  applicable  local  stationary,  homo¬ 
geneous  regimes. 


7.2  LOBD  Algorithms 

Locally  optimum  Bayes  detector  (LOBD)  algorithms  for  coherent  and 
incoherent  detection  are  found  under  the  above  detection  modes  (7.3a), 
(7.3b)  to  be  explicitly  [4]: 


A.  Coherent  Detection 


Sj’Wcoh  ’  BJ*-coh  -  1  j/<Vn><6m,n> 

m=l  n=l 


J  =  MN  (7.4) 


wi  th 


<0  >  =  <a  s^>; 

m,n  on  n 


=  37  109  "i<xlH0l 


x=x 


m,n 


<  =  Vn>' 

(7.4a) 


B.  Incoherent  Detection 


M  N 


^J^inc  BJ-inc  +  2!  n  ^m.nSn'  ,n'  +  Sn.n^mm' Ann '  ^  <'in,nBm' .n^ 


m,m  n,n 


where  now 


<6  9  ,  ,> 

m,n  m  ,n 


,2  ~  (mm1 ) 

3  ni  |  i  |  p  i  , 
o  n-n  in-n 


(7.5) 
(7.5 a) 


with 
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n-n 


a  a 
on  on 


(7.5b) 


mi 


(m-m;)  E  <;(.);(•')> 

«*1  _  ■"»  1  I  *■»  V*  1 


n-n 


and 


0  1 
m,n 


d£ 

dx 


x=x 


cf.  (7.4a). 


m,n 


In  the  above,  s^  =  s^(tn)  is  the  normalized  signal  waveform  at 
time  t  at  the  sensor,  with  <s^m^>  =  1;  a^  = 


on 


aon  =  ao(tn>  here 


(signal  field  uniform  over  the  array),  with  aQn  =  AQ(tn)//2^,  in  which 
4  =  I  =( I '  +  fi )  is  the  sum  of  the  intensities  (at  the  receiving  aperture)  of 
the  nongaussian  interference  I*  and  gaussian  noise  N.  (Here 


because  of  the  postulated  homogeneity  and  stationarity,  cf.  above.)  The 

normalized  data  samples  are  x  =  Xj_m  n//F  ,  with  the  desired  signal 

vector  S  =  {S^}  =  .  Here  w.  is  the  first-order  pdf  of  the 

total  interference:  because  of  independent  sampling  w,  =  n  w  (xm  ). 

u  m  ^  n  l  i  II  y  n 


7.3  Signal  Structures  for  Coherent  and  Incoherent  Reception 

Specifically,  we  note  also  that  in  realistic  situations  involving 
fading  and  doppler  uncertainties  the  desired  signal  structure  is 


t  =  nAt 
n 

(7.6) 

where  uj  =  random  doppler  shift,  e  =  e0  =  signal  epoch  which  is  a  priori 

known  in  coherent  detection;  <J>  =  possible  phase  modulation,  for  tne 

typical  narrow-band  signals  used  in  most  acoustic  environments.  The  net 

t  h 

path  delay,  At  ,  from  the  m— - element  to  some  selected  reference  point 
associated  with  the  array,  is  given  by 


(7.7) 

where  _Tq,  are  respectively  the  unit  vectors  in  the  direction  of  the 
(incoming)  signal  wavefront,  and  the  main  axis  of  the  beam  formed  by 


At 


m 


{lo 


r  /  c 
~m  o 


=  (v 


'W 
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the  equivalent  spatial  alignment  of  the  M  elements;^  is  the  (vector) 


J.  L- 

distance  from  the  element  in  the  array  to  the  reference  point.  Thus, 
when  the  beam  is  steered  to  the  signal  direction,  T  D  =  T  ,  cf.  Fig.  4.1,  and 

/  \  -  *0  K  '0 


At  =  0,  so  that  s„ 
m  n 


(m)  = 


=  sn:  the  sampled  signal  component  is  now  independent 


of  m  in  (7. A),  (7.4a)  above,  for  coherent  detection. 

For  the  incoherent  cases  (7.5),  (7.5a,b),  (where  eQ  is  now  random), 
we  have,  for  a  fixed  signal  wavefront  direction, 


Am')  =  e 
p I n-n  '  I  e 


■[Aci)cj(n-n 1  )At]*V2 


cos[cj  (n-n')At  +  k  (t  -1  D)*Ar  ,] 
L  o  o^o  ~-oR  — mm 


ko  ■  2*'xo  =  “o/co 


(7.8) 


with  Ar  ,  h  r  -r  1 ,  and  there  is  no  degradation  of  the  beam  pattern.  This 
~~~ mm  m  ~  m  /  . , 

signal  correlation  function,  P  j  n_n * |  »  is  clearly  maximized  (when  summed 
over  m,m',  in  (7.5),  for  example)  when  the  receiving  beam  is  pointed  at  the 
normal  (j  )  to  the  incoming  signal  wavefront,  e.g.,jQ  =  TqR.  Thus, 

Pi  n,,,  also.  Alternatively,  when  the  desired  signal  wavefront 


(mm1 ) 
n-n‘ 


is  perturbed,  arriving  from  various  different  directions,  beam  structure 
is  destroyed  and  P [n-n ' i  is  much  reduced>  with  a  consequent  serious  degrada¬ 
tion  of  performance  [cf.  (6.10b)  and  Sec.  8  ff.].  Doppler  smearing 
(A<jOj>k  =  0),  of  course,  always  reduces  the  wavefront  (and  waveform) 
coherence  of  the  desired  signal  [cf.  (7.6),  (7.8)],  similarly  degrading 
performance,  as  one  would  expcet  [cf.  Sec.  6.1  above]. 

Finally,  we  note  that  because  of  the  postulate  of  independent  sampling, 
space  and  time  processing  are  interchangeable  operations,  cf.  (7.4),  (7.5): 
we  can  process  over  the  array  (m=  1,...M)  at  any  given  instant  (tp),  or 
equivalently,  at  each  array  element  (m),  over  time  (n=  1,...N).  The  choice 
is  a  matter  of  technical  convenience.  In  addition,  we  observe  that  for 
coherent  detection,  (7.4),  (7.6),  (7.7),  single  beams  are  formed  (-s^  ;), 
while  in  incoherent  detection  product  beams  are  generated  (m,m‘)  in  (7.5). 


7.4  The  Bias  and  Asymptotic  Optimality  (AO) 

As  explained  above,  cf.  (7.3a,b),  the  bias  terms  for  are  found  here 

to  be 
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BJ-coh 


I  <a'!!)>2  <s'm)>2 


2  ~on 

m  i  n=i 


generally.  With  fading  and  doppler  smear  (7.9)  becomes 


'  ■»  n\ 


6J-coh  log  p  =  6J.coh  i  -l\  (l-'l)MNL(2,H1(NAta„d), 


(7.9a) 


2  2 

where  0<l-n^Cao>  /<a^><  1  is  a  fading  parameter  (no  fading:  n  =  0;  deep 
fadi ng:  n  =  1) ,  and 


H.(x)  =  /irG(x)/2x,  with  O(x)  =  (2//r)/  e  t  dt  =  erf  x,  (7.9b) 
1  0 

represents  the  effects  of  doppler  smear. 

Similarly,  for  the  incoherent  cases  one  gets 


B*  - 
J-inc 


J  ((L(4)-2L(2)2)4, 


m,m'  n,n' 


mm’ "nn 


,+2L(2)2,c(mn,^)|i;(m.'n')2 


n-n 


n-n 


.  <a(m)a^')>  <s(m)s(m,)> 
oo  n  n 


(7.10) 


which  reduces  for  rapid  (one-sided)  fading  and  negligible  doppler  to 


(7.10a) 


(For  no  or  slow  fading  one  sets  n  =  0  in  (7.10a).)  The  result  (7.10a) 

(41  (2r 

assumes  N»1  (M>1)  and  the  usual  nongaussian  noise  situation  Lv  '  =  2L'  , 

where  specifically  [4],  [16] 

(7.11a) 

(7.11b) 


Thus,  L^,  are  statistics  of  the  (total)  interference,  I,  including 
the  gaussian  noise  component. 


e  <£2>H  >  all  m,n,  =  J  (wj/w^w^dx,  >  0,  cf.  (7.4a); 

0  -oo 

e  /  (wj'/Wj^Wjdx  (>  0). 
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For  asymptotic  optimality,  ("AO"  as  or  J>>  1)  a  sufficient  con¬ 

dition- -always  fulfilled  for  small  (>  0)  input  signals,  e.g.,  1  »  a2  >  0-- 
is  that 

varHQ  9J  E  “  <9?0  5  (°oj)  =  "2^J  varH1  9J^  ^7-12^ 

As  noted  above,  this  permits  us  to  terminate  the  expansion  (7.1)  in  the 
form  (7.2),  for  large  J ( =MN ) ,  and  still  ensure  the  desired  small  prob¬ 
abilities  of  decision  error  in  the  decision  process  (7.3c).  [An  equivalent 

(sufficient)  condition  for  AO  is  that  gt,  as  J  -*■  ®,  is  gaussian,  with 
2  z 

variance  (a* j)  and  means  log  y  +  (o*j)  / 2  under  HQ,  respectively,  cf. 

Sec.  2B  of  [16].]  (In  most  earlier  work  the  proper  bias  was  omitted, 
cf.  Secs.  I I b ,  IV,  V A  of  [16],  so  that  the  resulting  detection  algorithms 
are  neither  locally  optimum  nor  AO.) 

7.5  Suboptimum  Algorithms 

Suboptimum  algorithms  are  handled  in  the  same  general  fashion. 

Instead  of  the  quite  complex  characteristics  (-8.  )  in  g^,  (7.4),  (7.5) 

one  can  use  much  simpler  forms:  for  the  simple  correlator  one  sets  8  n-*-x  , 

m  9 II  ill  5  n 

and  for  the  hard-limiter  correlator,  8  +  sgn  xm  „  (with  the  8'  „  term 

m,n  m,n  m,n 

omitted  in  (7.5)).  The  former  is  well-known  to  be  (threshold)  optimum 
generally  in  gauss  noise,  w.(x)  =  e"x  7  /  /2tT,  while  the  latter  is  threshold 
optimum  in  Laplace  noise,  w^(x)  =  (/2)  iexp(-/2| x| ) ,  <x  >  =  1;  <x>  =  0, 
cf.  Sec.  Ill  of  [11]. 

8,  Threshgl^De  tec  tor  ^Performance 

Because  of  the  large  samples  involved  (J  =  MN»1)  and  their  independence, 
the  test  statistic  g^  is  asymptotically  normally  distributed,  under  both 
Hn  and  H,  with  varu  g*  provided  the  input  signal  is  suitably  small  [cf. 

U  J.  ng  J 

(7.12)  and  II, B  of  l 16] ,  etc.].  This  permits  a  direct  calculation  of  per¬ 
formance.  For  the  important  detection  cases  here,  we  have  the  familiar 
Neyman-Pearson  (N.P.)  form  [15]  of  the  probability  of  correctly  detecting 
the  presence  of  a  desired  signal: 
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P0*}  =  f  {1  +0[  -  0-1(  1-2o4*})]}  .  (8.1) 

/2 

( *\ 

where  a p  '  is  the  (conditional)  false  alarm  probability  of  incorrectly 
deciding  the  signal  is  present  when  only  noise  occurs,  again  O(x)  =  erf  x, 
cf.  (7.9a),  and 


a 


(*) 

oJ 


For  suboptimum  detectors  we  drop  the  superscript  *. 
given  by 


a 


(*) 

F 


r~(*) 

oJ  ,  loq(K/u) 


2/2  /2  a 


oJ  J 


(8.2) 


Expl ici tly ,  a 


(*) 

F 


is 


(8.3) 


for  u,  cf.  Eq.  (7.3)  et  seq.  In  the  usual  applications,  y=l,  i.e., 
p =  q =  4-  (In  communication  applications,  when  the  decision  criterion  in 
the  so-called  Ideal  Observer  [15] ,  [16],  p =  1  also,  i.e.,  p  =  q  =  l2  and  one  has 
a  symmetric  "on-off"  channel  for  which  the  error  probability  (per  decision 
or  per  baud,  symbol,  etc.)  is  now 


[y  =  1;  K  =  1]. 


(8.4) 


8.1  Performance  Examples 

The  key  element  in  measuring  performance  is  clearly  a'/,  cf.  (8.1)- 

(*\0J 

(8.4).  We  list  below  a  variety  of  specific  forms  for  o',  . 


.(*) 


A.  Coherent  Detection:  (aoj)coh  =  2a2(l-n)M.NL^Hj(NAtAu^)  =  2a2  H*Qh, 


oJ 

i  ^  V  "  At  AUy  , 


(8.5) 


cf-  [4],  [15],  [16]. 


^  MN  .  (4)  ..  (2)2 


8.  Incoherent  Detection  [4]:  (o*,).  =  2ad  I-—-  L^4'+2Ll“' 

-  oJ  me  o  8  o  N 


22 

H  2ao  ninc’ 


(8.6) 


where  is  the  signal  structure  factor,  which  in  these  stationary, 
homogeneous  cases  becomes  explicitly 
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Qm  =  1  + 


1  N 
I  y  < 

N  l  . 


^  2  2  *  * 

N  4  .  ml n-n 1 1 ^ I n-n 1  I  ’  «^>o  «hdR’ 

n,n'  *  1  * 


cf .  (7.8)  et  seq. , 

(8.6a) 


which  is  maximized  when  the  beam  formed  is  directed  at  the  desired 

signal  source.  _ 

? 

Here  n*  ^  is  the  space-time-bandwidth  product,  and  the  aQ  factor  is 
defined  as  the  minimum  detectable  signal,  <(ao)>min‘ 

Evaluation  of  for  various  conditions  of  doppler  "smear"  and  fading 
gives  results  like  [4]: 


(  a*  )2 
v  oJ'inc 


2-siaed,  slow-fading; 
coherent  signal  (Aco^O) 


=  a2  (MNL^2 V2)2  = 


KPinc 


1-sided,  no, 
slow;  coh.  sig. 


(8.7a) 


Kj’lnc 


1-sided,  very  rapid; 
coh.  sig. 


af{MN(l-n)L(2)/2}2. 


(8.7b) 


In  all  cases  here  we  note  that  a* j ~  /f  :  increasing  the  "space-time-bandwidth 
product"  by  increasing  the  number  of  both  the  independent  spatial  and  temporal 
samples  (M,N)  clearly  improves  performance,  from  (8.1),  (8.4). 

It  is  usually  convenient  to  relate  the  probability  controls,  here 

(*\  (*\  .  2^* 

Pg  ,  to  the  processing  gain,  and  minimum  detectable  signal,  <a  >m|n, 

by  means  of  the  following  relations,  obtained  from  (8.5),  (8.6)  in  (8.7), 

(8.4),  viz.: 


,2><*>  . 
0  min-coh 


■<*)  ,2/n(*). 


<a‘>' . ’  ...  ■  (CJ  '  )‘/n 


coh  ’ 


<a2><*> 


(* 

=  r' 

0  min-inc  N 


k)  /M*T 

.P/y  inc  ’ 


(8.8) 


where  ^  p  =  0  1  ( 2pg  ^-1)  +  0  *(1-2  p^)  =  0qjV'/2  specifically.  From 

(8.8)  vs.  (8.5),  (8.6)  we  find  FI*  at  once. 

For  nongaussian  noise  fields  L^,  can  be  quite  large  0(10^,10^), 

-5 

typically,  for,  say.  Class  A  noise  with  A^  =  0.3,  F^  =  10  .  Similar  orders 
of  magnitude  are  noted  for  typical  Class  B  cases;  ([9];  see  Figs.  4-7  of  [16]). 
However,  for  gaussian  noise  =  1,  =  2. 

.(*> 


Figure  8.1  shows  the  probability  of  detection,  pi  ',  versus  the 

2  (*)  .  .  .  (*) 

pS'  ’  i.(  ,  +■  h  n^nroc  c  1  nn  n  a  ,  n  H'  ' 

o  min' 


minimum  detectable  signal,  oSi-',  with  processing  gain  FT  '  as  parameter, 
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<al  '>  . 
o  min 

Figure  8.1  Probabi 1 i ties  P i  ^  of  detection  (8.1)  for  both  coherent  and 
incoherent  reception,  with  processing  gain  as  parameter,  for 
false  alarm  probability  '  =  10'^,  10~4,  10"$. 


/  &  \  __  n  ^  n  r 

for  various  levels  of  false  alarm  probability:  ctp  '  =  10  ,  10  ,  10  . 

Clearly,  coherent  detection  is  better,  by  a  factor  2  in  dB,  than  incoherent 
detection  with  the  same  processing  gain,  cf.  (8.5)  vs.  (8.6).  For  example, 

-•5 

with  a  processing  gain  of  80  dB,  -  10  ,  and  a  required  pg  =  0.80,  we 

see  that  <a^in.coh  -  -33  dB,  while  <£*jn.inc  -  -66  dB. 


8.2  Remarks  on  Suboptimum  Performance 

The  performance  of  suboptimum  systems  is  readily  obtained  canonically 
(in  the  threshold  regime)  if  we  define  a  degradation  factor  q  .  by  the 


Here  n  is  thp  processing  gain  for  a  suboptimum  detector,  employed  against  the 
same  (not  necessarily  gaussian)  noise.  This  factor,  4^,  measures  by  how 
much  the  minimum  detectable  signal  for  the  suboptimum  system  must  be 
raised  vis-a-vis  that  of  the  optimum  to  achieve  the  same  performance, 

PQ  =  P*  (or  pD  =  p*),  viz.  oQj  =  a*j  in  (8.1),  so  that 


<atf*  min-coh  <aa>  min- 


coh^d-coh 


;  <aV*.  • 

a  mm-i 


=  <aT>  •  •  A a  •  ■ 

nc  0  mm-inc  d-inc 


(8.10) 


For  example,  for  a  coherent  detector  optimized  against  gauss  noise, 

=  1  in  (8.5),  so  that  4>d  =  (L^2;)~*«l  measures  the  degradation 

when  such  a  detector  is  used  in  strongly  nongauss  noise  instead  of  one 

-5 

properly  "matched"  to  this  noise.  Thus,  if  FI  =  10  ,  A.  =0.3,  we  have 

H  a 

-  50  dB,  Fig.  4,  [16],  so  that  the  optimum  detector  can  detect  a 
minimum  detectable  signal  50  dB  less  than  the  correlation  detector  in 
this  noise.  In  a  similar  way,  two  suboptimum  detectors  can  be  compared: 
4,(1  vs‘  2^  =  ri^1  Vri^  =  4)^/4>^2^.  The  quantity  is  also  seen  to  be 
equal  to  the  Asymptotic  Relative  Efficiency  (ARE,  9>0).  [See  Section 
IV,£  of  [16].  ] 


8.3  Discussion:  Threshold  Detection 

In  the  preceding  sections  we  have  sketched  some  of  the  main  results 
for  the  development  of  optimum  threshold  detection  in  nongaussian  and 
nonuniform  EMI  fields.  This  extends  earlier  work,  [55 ]- [57 ] ,  for  example, 
to  include  in  addition  to  the  nongaussian  field,  the  often  important 
situation  where  the  fields  received  by  an  antenna  (array)  are  not  uniform 
over  the  array  elements,  basically  because  of  the  size  of  the  antenna 
vis-a-vis  the  wavelengths  of  the  received  fields.  In  addition,  we  have 
outlined  briefly  various  elements  in  the  construction  of  statistical 
physical  models  of  the  acoustic  field,  including  the  needed  pdf's. 

Our  analysis  has  postulated  independent  (noise)  samples  in  space 
and  time.  The  latter  are  practically  achieveable  without  much  diffi¬ 
culty  by  slightly  increasing  the  sampling  period  (0(50%)).  Independent 
spatial  sampling,  however,  is  not  generally  possible  with  non-sparse 
arrays,  because  of  the  non-uniformity  of  the  noise-field  and  unequal 
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spacing  of  the  zeros  of  the  correlation  functions  of  this  field,  cf. 

(6.10),  (6.10a,b).  In  practice,  this  means  that  effectively  one  may 
have  1<M'  <  (<  )M  equivalent  "independent"  spatial  samples:  some  spatial 
processing  improvement  over  a  single  element  (or  single  beam,  M=l)  is 
obtainable,  but  not  the  full  M-factor  theoretically  achieveable  if 
independent  spatial  sampling  is  truly  obtained.  Nevertheless,  for  such 
noise  fields  as  (6.10)  (in  addition  to  improving  beam  resolution),  it 
can  be  worthwhile  to  use  relatively  large  arrays  (vs.  XQ),  as  M'  >  1, 
and  o*j,  -  /M'N' ,  cf.  (8.5) ,  (8.6) ,  is  larger  than  a*N  =  /N,  M  =  1. 

Of  course,  algorithms  like  (7.4),  (7.5)  are  not  optimum  when  the 
sampling  is  not  independent,  but  they  are  comparatively  close:  pre¬ 
liminary  estimates  for  typical  Class  A  interference  indicate  that  only 
0(2-3  dB)  improvement  is  theoretically  obtainable  if  all  correlations 
were  taken  into  account,  i.e.,  if  one  uses  the  J-joint  pdf  Wj(xJ^ 
rather  than  w^x^.)  f°r  the  basic  Pdf  used  in  the  general  expansions 
(7.2)-(7.3b)  above.  This  in  one  compelling  argument  for  developing 

"independent-sample"  algorithms;  another  is  that  there  are  no  tractable, 
t  h 

or  known  J~-order  models  (J  >  1,2)  for  the  underlying  pdf,  Wj . 

It  is  instructive  to  observe  that  beam-forming  is  automatically 
achieved  on  setting  the  relative  path  delays  Axm,  (7.7),  equal  to  zero, 
regardless  of  the  particular  geometrical  locations  of  the  array  elements. 
Thus,  we  can  have  arbitrary  distributions  of  elements,  but  still  recon¬ 
struct  the  desired  signal  wavefront  by  proper  delays  (i.e.,  maximize 
^(mm  )^  (7.8)).  Beams  formed  here  in  this  manner  are  adaptive,  weighted 

beams.  The  important  trade-offs  are  coherent  space-time  signal  wavefronts 
vs.  incoherent  (uncorrelated)  noise  field  samples,  thus  increasing  the 
effective  statistical  sample  under  vs.  H^,  i.e.,  lowering  the  (here) 
joint  probabi 1 ities  of  Type  I  and  II  errors,  e.g.  ,  m]n  (qa*+p6*),  •>  = 
decision  rule,  viz.  (7.3c). 

In  a  broad  sense  these  optimum  algorithms  seek  first,  thus,  to 
"match"  the  receiver  to  the  noise,  via  ^  n,  cf.  (7.4),  (7.5),  and  then  to 
"match"  the  desired  signal  to  these  processed  input  data,  2(xm  n).  It 
is  the  "tails"  of  the  noise  pdf  which  are  signficiant  vis-a-vis  gaussian 
noise:  the  nongaussian  interference  falls  off  much  more  slowly,  so  that 
the  transfer  characteristic,  i,  between  x-input  and  2-output  includes  a 
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limiting  or  clipping  feature,  as  well.  For  this  reason,  limiters  can 

be  good  approximations  to  the  actual  characteristic,  giving  performances 

not  many  dB  worse  than  the  theoretical  optimum.  Of  course,  actual 

differences  will  depend  on  the  particular  EMI  field  to  which  the  receiver 

is  adapted,  i.e.,  for  which  the  field  parameters  (ft^,  A,  r ' ,  a,  etc.) 

have  been  estimated.  Finally,  it  is  emphasized  that  threshold  detectors 

of  the  above  types  are  useful  and  give  noticeable  improvement  over  con- 

(*) 

ventional  systems  when  the  space-time-bandwidth  product,  H;  ,  is  large, 

t*)  J 

so  that  the  required  small  probabilities  of  error  '  can  be  obtained. 

(For  a  discussion  of  asymptotic  optimality,  AO,  bounds  on  etc., 

as  well  as  some  numerical  examples,  see  again  [15],  [16].) 

9,._  .Threshold  Signal  Estimation 

As  is  well-known,  in  the  Bayes  theory  [20],  [58]  estimation  is  the 
"twin"  of  detection:  the  former  aims  to  determine  the  particular  charac¬ 
teristics  of  a  desired  signal,  e.g.,  amplitude,  waveform,  phase,  frequency, 
etc.,  once  the  signal  has  been  detected,  while  the  latter  is  concerned 
with  the  basic  question  of  the  desired  signal’s  presence  or  absence.  In 
a  very  broad  sense  "estimation"  can  be  regarded  as  an  extension  of  de¬ 
tection:  both  subsume  appropriate  cost  functions  for  the  derivation  of 
optimal  algorithms  and  for  the  measures  and  comparisons  of  performance. 

In  fact,  in  a  more  general  Bayesian  sense  detection  and  estimation  are 
the  two,  coupled  component  elements  of  the  composite  process  of  signal 
extraction:  detection  and  estimation  joined  together  by  a  suitable  cost 
function,  reflecting  the  frequent  situation  where  it  is  not  completely 
certain,  for  the  estimation  process,  that  the  desired  signal  is  present, 
which,  in  turn,  results  in  biased  estimators  [59]. 

Here,  however,  we  shall  assume  that  the  desired  signal  is  known  a 
priori  to  be  present,  and  that  certain  waveform  and  parameter  features  of 
the  signal  are  to  be  estimated,  when  the  accompanying  acoustic  inter¬ 
ference  is  highly  nongaussian  (the  "classical"  situation  of  gaussian 
noise  is  a  special  case  of  our  general  model),  flow  the  estimators  are 
unbiased.  In  any  case,  we  may  expect  the  estimation  process,  particularly 
optimum  estimation  (in  the  3ayes  theory  [58],  [20]),  to  be  closely  related 
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to  the  detection  process,  since  both  are  derived  from  approximate  likeli¬ 
hood  functions  [60],  [61].  This  will  be  seen  explicitly  in  the  canonical 
threshold  theory  outlined  below  (cf.  Sec.  9.1). 

From  the  practical  underwater  acoustic  viewpoint,  signal  estimation 
is  a  common  communication  requirement:  signal  amplitude  (or  level), 
waveform,  frequency,  epoch,  modulation  are  each  important  elements  of 
the  reception  process,  whether  it  be  telecommunications,  radar,  tele¬ 
vision  (and  their  analogues  in  acoustical  and  optical  regimes).  Designing 
receivers  for  effective  threshold  performance  usually  ensures  better  (but 
not  necessarily  optimum)  performance  at  strong-signal  levels.  As  in  the 
detection  cases  [16],  an  optimal  theory  provides  limiting  estimation 
algorithms  and  performance,  which  are  models  to  be  approximated  in  prac¬ 
tice.  Also,  as  in  detection,  a  canonical  theory  is  possible  in  the  weak- 
signal  cases  (as  long  as  reasonably  large  time-bandwith  products  are 
permitted,  of  course:  the  desired  signal  must  be  "extractable"  under  the 
observational  constraints).  Such  a  canonical  theory  also  provides  standards 
against  which  practical,  suboptimum  (i.e.,  approximate)  algorithms  can  be 
evaluated  and  compared  [15],  [16]. 

As  in  threshold  detection  theory  [16]  a  fundamental  problem  now  in 
developing  effective  threshold  estimation  algorithms  is  to  obtain  ex¬ 
pressions  of  limited  complexity,  which  retain  their  optimum  nature  when 
sample  size  (J)  becomes  increasingly  large.  The  latter  is  the  case,  of 
course,  when  the  signal  is  weak,  since  large  effective  data  samples  are 
required  for  small  expected  errors  in  the  resulting  estimates--anal ogous 
to  small  probabilities  of  decision  error  in  optimum  (threshold)  detection. 
However,  just  as  in  detection  [16],  without  a  suitable  “bias"  term  in  the 
extraction  algorithm  this  al gori thm  demands progressi vely  more  terms  in 
its  approximative  form.  This  rapidly  defeats  the  key  requirement  of 
processing  feasibility,  particularly  for  signals  and  interference  of 
practical  use,  and  destroys  analytic  tractability  as  well. 

Accordingly,  we  parallel  our  preceding  summary  review  of  (binary)  space- 
time  threshold  detection  theory  (Secs.  7,  8),  and  present  now  a  concise 
summary  of  optimum  and  suboptimum  threshold  signal  estimation,  for  the 
three  generic  modes  of  reception  noted  earlier  for  detection:  (1)  coherent 
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estimation,  (2)  incoherent  estimation,  and  (3)  composite  estimation,  where 
a  linear  combination  of  the  coherent  and  incoherent  algorithms  is  employed. 
[Because  of  spatial  limitations  we  shall  confine  our  explicit  illustrations 
to  the  simplest  situation  of  amplitude  (or  scale)  estimation,  where  signal 
waveform  is  otherwise  known  at  the  receiver.  However,  the  effects  of  fad¬ 
ing,  doppler,  and  other  propagation  conditions  can  be  included,  as  in  the 
detection  examples  above,  Secs.  7,  8.  The  general  formalism  outlined  here, 
nevertheless,  is  quite  capable  of  providing  explicit  algorithms  in  the  more 
complex  situations  .involving  waveform,  frequency,  phase  estimation,  etc.] 

As  before,  the  interference  is  additive,  generally  nongaussian,  and  is 
either  Class  A  or  Class  B  interference  (i.e.,  respectively  coherently  or 
incoherently  received  in  the  (linear)  front-end  stages  of  a  typical  narrow 
band  receiver,  accompanied  by  (additive)  external  and  internal  gaussian 
noise  [7],  [15],  [16],  [18].  These  interference  models  account  for  most 
of  the  practical  cases. 

9.1  Canonical  Optimum  Threshold  Estimator  Structures;  AQ/L03E  Forms 

A  canonical  optimum  estimator  is  represented  by  y*  =  y*(9|X),  where 
0.  =  (0  ,...,0M)  is  the  parameter  set  to  be  estimated  on  the  basis  of  j(; 
the  received  data,  a  refers  to  the  pdf  of  Q  (see  Appendix  A. 3).  As  noted 
at  the  beginning  of  Sec.  9,  the  critical  problem  in  developing  explicit 
locally  optimum  threshold  estimators  (LOBE's)  from  the  general  (all  signal 
level)  forms  (A. 3-6),  (A. 3-8)  is  to  terminate  the  desired  series  approxi¬ 
mation  properly,  i.e.,  without  an  excessive  number  of  terms  and  in  such  a 
way  that  these  estimators  are  asymptotically  optimum  (AO),  as  well.  This 
situation  precisely  parallels  that  for  optimum  threshold  detection,  already 
treated  by  the  author  in  [15],  [20],  In  fact,  the  approach  here  is  to  con¬ 
vert  the  various  likelihood  functions  oFn,  oWn  above,  from  which  the  LOBE's 
are  determined,  cf.  (A. 3-4),  (A. 3-8),  into  likelihood  ratios,  equivalent  now 
to  a  suitable  statistical  test  among  two  distinguishable  hypothesis  (H^,  H^), 
etc.  Then  leCam's  results  [60]  as  applied  by  Levin  [61]  are  used  to  obtain 
at  least  a  sufficient  condition  of  AO  for  the  "LOBD"  form  from  which  the 
resultinf  AO  LOBE  is  then  constructed,  according  to  the  error  function 
(QCF,  SCF,  etc.)  chosen. 

A.  Equivalent  LOBE's  for  LOBD's 

We  consider  first  the  SCF  case,  cf.  (A. 3-7),  with  the  optimizing  rule 
(A. 3-3)  for  generating  the  UMLE's,  Y*.  Since  the  (linear)  operation  of 
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differentiation  in  (A. 3-8)  is  unaffected  if  we  subtract  log  Fj(XjO),  we 
see  that  (A. 3-8)  becomes 


Hog  a(6s)  +  t(10)(X|8S)SCFi 


nr  SCF  Q,=yt 
m  'm 


=  0,  m  =  (9.1) 


where  0  =  ( 0~ ,9, 1 ) ,  in  which  0~  is  the  particular  parameter  to  be  esti 
mated,  with 

0  (10) 

“  1  n  dr  rvlc/a.  a'Us  /  r 


‘j!sCF  =  ’09i<Fj{X|S(6fi.9')}>9,/  Fj(X|0)] 


(9.1a) 


now  the  (logarithms  of  the)  generalized  likelihood  ratio  representing  the 
test  states  H^SeN)  vs.  HQ(N),  which  is  the  familiar  "on-off"  or  "signal 
and  noise"  vs.  "noise  alone"  situation  in  detection  theory. 

In  a  similar  way  we  see  that  the  rule  (A. 3-6)  for  obtaining  the 
Bayes  estimators  for  the  QCF  (A. 3-5)  can  be  equivalently  expressed  in 
terms  of  another  likelihood  ratio  by 

n(2l)  t  V  \a^\ 


<\1QCF  =  {/%-aSa(eA)e  ° 


m'QCF 


d0~}  , 


(9.2) 


where  now 


*J?0CF  5  1<J9  /<Fj(X[S»e 


(9.2a) 


This  likelihood  ratio,  represents  the  binary-signal  detection 
situation  where  we  test  H^S-ffiN)  vs.  (Hj(S®N)»  i.e.,  =  JS^,  that 
portion  of  JS  containing  0-  only,  and  ^  =  S,  cf.  (9.2a).  Thus  we 


-(10) 


note  that  y*  cr[-  is  a  linear  functional  of  ,  ' ,  cf.  (9.1),  whereas  the 

Ul*oLr  J  ( 0]  \ 

{ Y*_qqf^ are  (monotonic)  nonlinear  functionals  of  ?.  j  ,  cf.  (9.2). 

Having  cast  these  estimators  in  likelihood  ratio  form,  we  can  now 
proceed  directly  as  in  the  previously  developed  optimum  threshold  detec¬ 
tion  theory  ([13],  [16],  Appendix  A. 3;  [15],  Appendix)  to  obtain  the 
needed  approximate  forms  for  ij  =  9j(x_>9^)  +  tj(x_,e~)  =  g*j,  where  it  is 
shown  that  as  J  -*  <»,  t  j  ■*  0  prob.  1  on  vs.  Hq  (or  vs.  H^),  i.e., 
lj  =  is  AO  as  J  -*■*>,  cf.  (7.12)  et  seq.  Because  of  the  functional 
relations  (9.1),  (9.2),  ,*  is  likewise  then  AO,  as  well  as  LOBE,  i.e., 
is  AOBE  (asymptotically  Bayes  estimator).  Now  specifically  (9.1),  (9.2) 
become 
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SCF : UMLE 


3log(W 


3gJ-SCF 


m  'm 


(9.3) 


{y *}  -  e^letfe 


9j-qcf 


(9.4) 


where  the  g*j  are  given  explicitly  in  the  next  section.  A  sufficient  con¬ 
dition  that  gj  is  AO  is  that  gj'  be  asymptotically  normal,  with  means 

?[varo9J-SCF’varl9J-QCF^  "=  H0’  Hl;  or  Hl»  H2 ’  and  variances  varogJ-SCF’ 
varlgJ-QCF’  resPecti vely.  [Note,  also,  that  because  we  have  assumed  the 

desired  signal  is  a  priori  known  to  be  [55]  present,  logy  =  log(p/q)  of 
the  detection  forms,  cf.  (2.4)  of  [20],  are  omitted  here  in  the  forma¬ 
tion  of  above.] 

B.  AO  Forms  of  the  g^ 

From  Sec.  3.2  of  [15],  and  from  [20]  suitably  adapted,  we  have  spe¬ 
cifically  the  needed  structures  of  g^  required  in  (9.3),  (9.4)  for  de¬ 
termining  the  associated  AO  LOBE'S.  Specializing  to  the  practically 
useful  cases  of  independent  data  (i.e.,  noise)  samples  and  locally  sta¬ 
tionary  noise  processes,  which  we  shall  henceforth  assume  in  the  remainder 
of  our  treatment,  we  have  once  more,  cf.  Sec.  7.2 

( 1 )  Co  he  rent  Reception  (SCF) 


9J-coh 


ToJ-coh  "  ^<aojSj>9' 


(2)  Incoherent  Reception  (SCF 


2  oJ-mc 


l  (*iV'Lij)<aoiaojsisj>n. 

J  v 


(3) 


Composite  Reception  (SCF) 


(9.7) 


where,  as  in  Sec.  7  et  seq.,  =  aQ ( t^ ),  etc.  is  a  (normalized)  signal 
amplitude,  s.  =  s(t.  ,0,8')  is  a  normalized  signal  waveform,  such  that 

2  '  i  ^ 

<s  >  =  1,  and  £.  =  d  log  w- (x^)/dx^,  2/  =  di^/dx^,  where  w^x-)  =  first 

order  pdf  of  the  interference  (nongauss  +  gauss)  above,  in  the  usual  way. 

Specifically,  here  and  henceforth  in  Sec.  9,  a  ,•  =  a^ 

oj  on 


(m) 

s .  =  s'  ,  since 


j  =  m,n,  cf.  Sec.  7.2  et  seq.  In  this  way  we  include  the  spatial,  as  well 
as  temporal  description  of  the  sampled  waves.  Thus  i  =  (m,n);  j  =  (m‘,n'), 
when  used  in  the  incoherent  cases.  „ 

The  proper  bias  terms  are  given  by  B*  =  -cry  / 2,  cf.  (7.12),  or 

(7.9),  (7.10),  with  (7.11),  etc.  [We  must  be  careful  here  and  henceforth 
to  distinguish  between  the  various  averaged  signal  parameters  (01)  and 
the  unaveraged  ones  {e~},  which  are  to  be  estimated.] 

Similarly,  for  the  QCF  the  associated  g^  needed  in  (9.4)  are 

( 1 )  Coherent  Reception  (QCF) 


(9.9) 


(2)  Incoherent  Reception  (QCF) 


(9.10) 
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(9.11) 


B 


(21)* 

J-coh 


.(21)* 

:J-coh 


(2)  J 


j  [<aojSJ>0' 


<a  .s.>?] 
oj  j  0 J 


B;2i>*  - 
J-inc 


.(21)* 

7J-inc 


=-  |  I  {(L(4)-2L(2)2)6i  .  +  2L(2)2}ARir  (9.12) 


ij 


where  AR^  =  ^a^s^j,  -  <aoiaojsiSj>|  (^(Ap)^). 

C.  Suboptimum  Estimators 

We  parallel  the  treatment  in  [15],  [20],  where  g^  is  replaced  by 
the  suboptimum  gj,  now  appropriately  adapted  to  the  estimation  structures 
(i.e.,  hypothesis  test)  required  by  the  SCF,  QCF,  etc.  We  illustrate 
the  procedure  with  the  case  of  SCF  and  the  use  of  simple  (cross-  and 
auto-)  correlation  receivers,  well-known  to  be  threshold  optimum  in  gauss 
interference.  Now  (9.5)- (9-7)  reduce  to  this  case  directly  on  setting 


(2) 


(4) 


1 ,  L^['ss  =  2.  Similarly,  for  the 


li  *  -xi  <•••  lf  =  -n.  a"d  •••  Lgauss 

QCF  cases  (9.8),  (9.9),  we  make  the  same  substitutions,  in  both  the  bias 
and  data  portions  of  g^p- 

In  the  important  cases  of  clipper-correlators,  which  are  optimum  in 
Laplace  noise  and  which  are  well  known  to  be  effective  against  non- 
gaussian  noise,  cf.  [15],  we  replace  by  sgn  (and  omit  Ji! ) ,  with 
appropriate  modifications  of  the  bias  terms,  cf.  Sec.  3.3  of  [15]  for 
details.  The  basic  pdf's  for  which  these  various  correlation  detector- 
forms  (and  hence  estimator  forms  via  (9.1),  (9.3)) are  optimum,  are  re¬ 
spectively  given  by  w^(x)  in  Sec.  7.5  above. 

Finally,  it  is  important  to  note  that  these  LOBE's  [obtained  from 

2  2 

(9.1),  (9.3)]  are  only  practically  AO  and  LOB  when  9  (=<a;>)  is  sufficiently 

2  0 

small,  as  in  detection,  i.e.,  there  exists  an  <a  (<<  1),  with 

<aS<<a:>m!1  .  such  that  for  a£  >  <a:>  ,  the  threshold  optimal  character 
o  o  max  oo  max 

of  g^  breaks  down  and  becomes  suboptimum.  At  larger  input  signals  g^ 

{-*  g,  now)  and  hence  y* etc.,  may  or  may  not  be  monotonical  ly  better, 
in  absolute  terms,  than  g^j,  a2<a|  Usually,  for  sufficiently  strong 

signals,  g^  (now  gj)  is  absolutely  better,  unless  the  information-bearing 
portion  of  the  signal  is  destroyed  by  the  algorithm  itself,  e.g.,  clipping 
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destroys  waveform  detection/estimation,  but  not  phase  parameter  extraction, 
for  example. 


9.2  An  Example:  Amplitude  Estimation 
with  the  SCF:  Coherent  Reception 

The  simplest  useful  example  is  that  of  estimating  the  scale  (or  ampli¬ 
tude)  ~aQ  of  an  otherwise  fully  known  signal  waveform  received  in  (gen¬ 
erally)  nongaussian  interference.  This  is  also  the  problem  of  estimating 

~  ~  2 

signal  intensity  i*  =  a*  . 

Here  we  have  0~  =  0.  =  a  ,  with  0'  all  other  relevant  signal  parameters. 

With  coherent  reception,  signal  epoch  e  is  precisely  known,  e.g.,  e  =  e  , 

and  we  select  <sj>e  =  sj_max  =  with  independent  noise  samples. 

Thus,  we  write  <s.>Q(  =  s  =/2  here.  Physically,  our  present 
i  max 

example  can  represent  slow  fading,  whose  changes  are  negligible  over  the 
data  acquisition  period. 

Accordingly,  (9.5)  reduces  directly  to 

(gJ-cohW  =  'L^aoJ"ao/2J  (9‘12) 

Applying  this  to  (9.3),  with  0m  =  0^  =  aQ,  gives 


°l(ao>' 

Ol(ao> 


L<2)aoJ  -  V2  l  e(xi> 

<3 


=  0, 


a  =a* 
o  o 


(9.13) 


whose  solution  yields  the  desired  estimator  a*(x). 

To  proceed  further,  we  need  to  know  the  a  priori  pdf  of  9j(=aQ), 
e.g.,  '-jUq)-  bet  us  assume,  then,  that  J,(aQ)  =  l/AaQ,  a  uniform  pdf 
on  (0---'  aQ  <  aQ+),  and  zero  elsewhere.  Consequently,  a|(a0)  = 

0<ao<Aao  (the  contribution  of  the  ^-functions  at  0,  £aQ  provide  no 
meaningful  solutions).  Solving  (9.13)  gives 

a2<vl  (9.14) 


r 

a;<i> 

coh 

=  -  1  2-,/J^2L(2) 

SCF 

j 

uni  form 
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The  optimum  threshold  estimate  of  signal  intensity  is  accordingly 


(9.15) 


Here  §*,  and  •••  §*  are  UMLE  for  the  assumed  uniform  pdf  of  a  .  and  simul- 
o  o  o  _ 

taneously  CMLE  for  any  (meaningful)  pdf  of  a  .  Moreover,  <aT>,  =  a  . 

__  _  /  A  \  0  0  10 

because  <S>^  (=  =  -/2  aQ  L'^  ,  and  we  can  show  that 


var.a*  =  a^(l  +  0(\))  -  a  ^  =  var  a„ 
1  o  ov  ' '  o  o 


J»1  ’ 


(9.16) 


from  <£.£ .>.  =  2a  L 
i  J  1  o 


(i  f  j).  Thus,  a*  is 

_  m  _  I  * 

unconditionally  unbiased,  e.g.,  <a*>^  =  a  (  =  AaQ/2)  and  a*-*aQ,  ^rom 

(9.16). 


o 

The  "smallness"  condition  on  aQ  is  (from  Sec.  6.4,  Eq.  (6.71)  of  [15]) 


.  l(2)/[l(2.2)  .  a(2)‘]  ,  x#.  l(2.2)  , 

(9.17) 


var  V 
o 


2  2 

Here,  from  a, (a  )  =  1/A an  we  have  a„  =  (Aa  )  /3  in  (9.16)  to  establish 

1'  o'  o  o  _  o  ? 

upper  bound  (a£)  («  l)  for  which  a£  ^(af)  „  and  a*(x)  is  then  LOBE 

o  max  '  o  o'max  o  ~ 


an 
and  AO 


As  a  numerical  example,  consider  a  Class  A  noise  for  which  =  10 


-4 


A.  =  0.5,  so  that  x.  =  1.7*10~4  (from  Fig.  5.6  of  [15]):  •**  x  dB  =  -  37.7  dB, 

n  O  a  q  0 

with  the  result  that  a£«  =  -  38  dB  or  (Aa  )  << -37.7  +  4.8  =  -32.9  dB. 

V  V 

Typically,  we  might  choose  (a„)m,„  =  -43  or  -48  dB  here.  This  means  that 

O  flla  X  2 

sample  size  J  should  be  large  enough  to  ensure  (9.16),  i.e.,  J  =  0(1(J  or 
more)  appears  sufficient. 


9.3  Discussion 

First,  let  us  note  that  our  illustrative  example  is  readily  extended 

to  the  estimation  of  waveform  itself  JL  =  {S-}:  9~->-S~,  and  we  get  a  series 

W  J  m  m  3 

of  relations  like  (9.13),  if  o(£)  =  ^5^(9-);  otherwise  one  has  M-coupled 
equations  to  be  solved  for  0^  =  S~.  Also,  if  the  QCF  is  used,  the  structure 
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of  the  optimal  estimator  is  generally  much  more  complex,  due  to  the  non¬ 
linear  functional  relation  between  and  Yj>  cf.  (9.2).  If  the  parameters 
to  be  estimated  appear  functionally  in  the  waveform  S(t,£)  then  one  ob¬ 
tains  functional  solutions  for  the  estimators. 

The  important  general  result  here  (for  all  signal  levels)  is  the 
explicit  functional  relation  between  the  (optimum)  detection  algorithms 
and  the  resultant  optimum  estimators,  cf.  (9.1),  (9.2):  detector  structure 
provides  the  initial  analytic  relation,  from  which  the  estimator  is  the 
appropriate  functional  relation,  the  form  of  which  depends  on  the  choice 
of  cost  function.  Thus,  to  carry  out  the  derivation  of  the  desired  es¬ 
timator,  we  must  start  with  the  appropriate  detector  algorithm  [16].  The 

evaluation  of  estimator  performance  is  provided  by  using  the  estimator, 

(*) 

Yj  ,  in  the  expressions  for  the  average  error  (or  risk),  cf.  (A.3-2). 

Finally,  one  important  extension  of  the  analysis  outlined  above  is 
to  the  frequently  occurring  situation  (mentioned  at  the  beginning  of 
Sec.  9)  necessitating  joint  detection  and  estimation,  when  it  is  not 
known  precisely  that  the  desired  signal  is  present  [59]. 


10.  Concluding  Remarks 


In  the  preceding  section  we  have  outlined  a  comprehensive  methodology 
for  handling  the  general  problems  of  acoustic  signal  transmission  and 
reception  in  and  through  complex  underwater  media.  Our  approach  has  in¬ 
cluded  the  following  important  channel  features:  (1)  inhomogeneous  volumes 
and  interfaces,  both  deterministic  and  random;  (2)  spatial  as  well  as  tem¬ 
poral  effects;  (3)  the  inherently  nongaussian  character  of  major  components 
of  the  random  fields  arising  both  from  ambient  sources  and  produced  by  a 
variety  of  scattering  mechanisms;  (4)  space-time  processing  by  the  dis¬ 
tributed  arrays  which  couple  source  and  receiver  to  the  medium;  and  (5) 
threshold  detection  and  estimation  of  signals  in  these  generally  nongaussian 
environments . 

In  addition  to  the  above  topics  we  have  touched  upon  a  variety  of 
specific  methods  and  techniques  for  channel  modeling,  including  operator 
formalisms  and  diagrammatic  methods.  A  number  of  new  results  are  also 
discussed:  poisson  field  models  and  statistics,  canonical  threshold 
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estimators  and  estimates,  for  example,  as  well  as  measures  of  threshold 
detector  performance  which  specifically  incorporate  the  rfiles  of  the 
receiving  array,  signal  structure  (with  doppler  smearing  and  fading),  and 
mixed  gaussian  and  nongaussian  noise  components  in  the  interference. 
Spatial  processing  gains  are  seen  to  be  achieveable  in  the  usually 
encountered  nonuniform  noise  fields  by  suitable  spatial  sampling,  for 
both  detection  and  estimation,  in  addition  to  the  capability  of  beam 
formation.  Suboptimum,  as  well  as  optimum  processing  algorithms  are 
noted. 

The  preceding  sections  are  intended  to  be  a  general  guide  to  the 
formulation  and  treatment  of  specific  propagation  and  signal  processing 
problems,  wherein  attention  is  called  to  the  many  physical  features  which 
may  have  to  be  taken  into  account  in  individual  problems.  The  accompany¬ 
ing  references  should  provide  the  in-depth  analysis  and/or  specific 
results  needed  at  each  stage  of  the  problem  in  question.  As  stated  at 
the  beginning,  our  general  aim  here  has  been  to  provide  an  overview, 
keyed  to  the  realities  of  specific  applications,  and  at  the  same  time 
to  exhibit  a  unifying  and  interdisciplinary  approach,  within  which  the 
many  special  modeling  and  processing  tasks  of  underwater  acoustic  inter¬ 
ests  are  embedded.' 
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Appendix  A.l  Diagram  Representations  and  Solutions  [2] 


Diagram  representations  and  methods  are  very  useful  and  compact  forms 
for  describing  field  interactions  as  is  well  known  [23],  [41 ].  For  some 
of  our  subsequent  results  we  shall  include  diagram  equivalents,  along  with 
the  corresponding  analytic  expressions,  defining  the  symbol  "vocabulary" 
in  a  consistent  way.  For  our  ensemble  results  [(3.4)-(3.9)],  we  write 
the  equivalent  ensemble  diagrams: 


FD: 


r~— 1  =  177777773  - - »-*—  ol - 1:  (Eq.  (3.4); 

(a)  aH  MJQ}  ,{a} 


(A. la) 


(A. lb) 


Fos:r  rr — 3 

Eq.  (3.8) 


— - -  17777773 

1— *-o  H 


\mum-- -mjg 


H 


(A.lc) 


PTSS : 


CZZJ  = 


+  — ^  *r~  O- 


.J  C22ZZZ23 ;  Eq. 
H 


(3.9), 
(A. Id) 


where  we  define  the  symbols  as  we  go  along.  Thus,  ■  ->  (eM^)  denotes 
the  "feedforward"  operator,  while  « — o(e{Q})  denotes  the  (ensemble) 
"feedback"  or  scattering  operator.  The  ensemble  Feynman  diagram  (FD) 
equivalent  of  the  FOR  (A. la)  is  just  (A. lb),  with  (A.lc)  giving  the 
corresponding  ensemble  FOS,  and  (A.lc)  the  PTSS,  obtained  by  iteration  of 
(A.  lb)  or  from  (3.9)  directly. 

A. 1-1:  Stocnast ic  Solutions 

As  we  have  noted  above  [cf.  Eq.  (2.8)  et  seq.]  the  solutions  to 
the  Langevin  equation  (2.3)  are  the  various  statistics  of  the  field  t, 
e.g.,  here  the  moments  <c>,  <  ^  t£> ,  etc. 


A-l 


Let  us  consider  first  the  mean  field  <ct>.  Taking  the  average  of 
(3.9),  for  example,  gives  directly 

<>>  =  aH  +  l  <T|(n)>an  ,  (A. 2) 

M  n=l 

i,  1_ 

which  shows  that  all  (n—  -)  order  moments  of  Q  are  required.  The 
diagram  equivalent  of  (A. 2)  is,  from  (A. Id), 

[<=»=  [i+  r-&\+  +  ...]  I/Z////X  ;  with 

)  <o>  M<Q> 1  m1<Q1m2Q2>  aH 

•  =  <Q>;  \  =  <Q1...Q>.  (A. 2a) 

(n)  1  n 

There  is  a  second  approach  which  leads  to  an  equivalent  determin¬ 
istic  formulation,  which  is  usually  more  effective  for  approximations. 
Defining  a  deterministic  operator  qjd^  by 

q[d)<o>  E  <Qa>  ,  (A. 3) 

we  may  again  average  (3.4)  and  use  (A. 3),  to  get 

<o>  =  aH  +  MQ[d)<a>  =  (1-Tj (d)  )‘\  =  +  l  (T(  (d } ) (n ^ ,  (A. 4) 

or,  in  diagram  form. 


4=^  =  G777ZZZ1+ 

H 


—  7777ZZ2 

!—*«-•  H 


=  [1+  I  ( - M-«)(n)]EZZZZ3  ) 

1  H 

(A. 4a) 

respectively. 

Equation  (A. 4)  is  a  form  of  Dyson's  equation  (DE),  where  now  the 

y  /  ,  v 

equivalent  deterministic  scattering  operator,  Q'  or  (EDSO),  is  the 
analogue  of  the  "mass-operator"  (Q)  in  quantum  electrodynamics  (cf.  *160, 


A-2 


[24];  [23]).  The  EDSO  is  always  an  integral ,  or  global  operator,  over 
all  pertinent  space  and  time,  and  has  the  form 

ojd)  =  :  V<Q>  i  Am=  <QMQ(m)>-  Vfim.j <(MQ)(j)>, 

(A. 5) 

which  at  once  permits  a  hierarchy  of  approximate  forms,  by  stopping  with 
term  (m>l).  The  diagram  equivalent  of  (A. 5)  is 


i(d)  = 


[•]m=0+[4*  -  A-*]m=1+[4*£*  -  4W-  +  ^-4^4=2 


(A. 5a) 


Similar  concepts  may  be  applied  to  higher-order  moments  of  the  field. 
In  fact,  we  may  write  the  following  FOS  for  the  second-order,  second 
moment  of  the  field: 

<0^01 2>  =  (i-r]^))*1(i-r)}d)r)^d))<a1>  <^2> i  r]{d)  = 

(A. 6) 

where  the  EDSO,  Q^d^,  is  defined  by 


^12  <aIa2>  S  <Q1Q2aia2>' 


(A. 6a) 


cf.  Q|a  ,  (A. 3)  above.  Equation  (A. 6)  is  a  second-order  Dyson  equation, 
analogous  to  a  form  of  the  Bethe-Salpeter  equation  in  quantum  electro¬ 
dynamics  ([21];  1160,  [24]).  Analogous  to  (A.4a),  we  have 


=  EZZZZ22  + 

12  12 


<a^ot2> 


mm 


12  12  12 


r  i 


;  Eq.  (A. 6a); 


1  2 


,  etc .  ,  (A. 7a ) 
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with  obvious  extensions  of  the  second  and  third  diagrams  in  (A. 4a).  The 
EDSO,  (A. 6a),  is  the  extension  of  (A. 5a): 


d  / \  d  d 

O'd  J  =  < — •  =  r  4  b  _  •  •  +<  »  <  -•  -)  + 

q12  12  L  1  2  1  2  12  Jm=0 


(A. 7b) 


Higher  order  Oyson  forms  of  the  FOS,  like  (A. 6),  can  be  constructed 
for<a....a  >.  In  addition,  various  approximations  can  be  instituted,  in 

v\  /  _j  \  h* 

Cra;  (series  modification)  where  the  series  are  infinite  (for  some  form  of 
strong  scattering),  or  by  truncation,  for  "weak"  scattering  (lst-Born 
approximation,  for  example  [cf.  (5.1)]),  where  one  stops  at  the  dotted  line 
in  (A. Id),  (A. 2a),  with  n =  1  in  (A. 4a),  etc.  Unfortunately,  space  does 
not  permit  us  to  pursue  these  topics  further  here. 


A. 1-2.  Received  Waveforms 

The  expression  for  the  received  waveforms,  after  sensing  by  the 
receiving  array  (R)  may  be  similarly  diagrammed.  We  start  with  (A.l), 
etc.  applied  to  (4.1)  and  (4.2).  Thus,  we  have 


d 


(A. 8a)  and,  in  PTSS  form,  cf.  (A. Id), 
Eg.  (4.2): 


etc.  cf.  (A. 4), 


11  22 


11 


*-•]  — y 

22  1 


11 


+  — 

2  ;  11  22  22 

22  22  22  11  11 


22  11  11 


12 

(A. 8b) 


The  diagram  of  the  mean  field  is  given  at  once  by  (A. 2a),  or  (A. 4a), 
on  multiplying  the  receiver  operator  R  since  <X>  =  R<:<>.  The  diagram  for 
(4.5),  in  PTSS  form  for  <X^X^>  is,  similarly  by  direct  expansion  and 
averaging , 
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Eq.  (4.5): 


<X1X2>  =  +{— >(— : J*  +  j  +  +-+(-*+\  +  -  •  •  )rr  +  - 


11  22  1  1  2  2  2  1  1 


HEzzza 

H,„ 


(A. 8c) 


from  which  the  various  higher-order  moments  of  Q  (at  points  1,2)  appear 
clearly. 


Appendix  A. 2  Canonical  Scattering  Operator  formal i sms : 

Basic  Decqmposi_tion  Principle 

For  this  it  is  necessary  now  to  use  a  canonical  formulation  of  the 

A 

scattering  operator,  Q,  and  in  particular,  the  author's  concept  of  a 
Basic  Decomposition  Principle  (BDP),  which  permits  both  a  statistical  and 
physical  quantification  of  this  operator  in  terms  of  differential  scatter¬ 
ing  elements  (dse's)  which,  in  turn,  can  be  identified  with  specific 
physical  mechanisms.  Again,  space  limitations  allow  only  the  briefest 
summary  of  results  below. 

We  write  the  following  canonical  expressions 

oo 

Q-^  QM(X,t|r_' ,t' )  =/  dt‘  /  dN(Zir_'  ,t'  ;t)r1hM(t-t'  , tjjr '  ,r,. . . ) 


(A. 2-1) 


where 


dNf1  =  the  number  of  available  dse's  in  the  domain  d.l^;  /dN(  )  is  a 
"counting"  functional. 

h^  =  weighting  function,  rt  'esenting  the  reradiation  response  of 
a  dse,  now  postulated  to  be  a  linear  filter; 
d  ,,  =  the  local  interaction  operator,  determined  by  the  Langevin  equa¬ 


tion  (and  boundaries),  relating  the  dse  and  incident  field. 


(A. 2-la; 


Here  Zn  is  the  space  of  all  radiation  "events,"  e.g.,  scatteri-ng,  and 
Z  ( =Zpx ^  the  total  random  variable  and  parameter  space  (£);  (r,r“) 
are  vector's  in  the  coordinate  system  of  the  scattering  elements,  In 

most  cases  we  use  the  basic  point-scatter  model  (BPSM)  for  h^: 

hM(f  ,t|jr' ,. ..)  =  Yo(_r'  ,t-t'  )M5(t-AtdS);  Atd$  =  f  (t 1  ,t;r_' ,0  ,  (A. 2-2) 

where  Atrf^  =  t^-a  is  the  displacement  (in  time  units)  due  to  doppler 
of  the  dse,  away  from  some  fixeJ  position  in  and  yq  (>  0)  is  a  "cross 
section"  which  depends  on  the  nature  of  the  medium  (or  boundaries). 

The  inhomogeneous  character  of  the  medium  (and  boundaries)  is  of 
two  types.  A,  B.  We  have,  with  defining  conditions: 


A.  The  Inhomogeneum:  (0/)A  |7u|+  fr^oU),  c0  =  x0f0’ 

u  Lo 


(A.2-3a) 


B.  The  Discontinuum:  Ajvy|+|-  |^->>o(l),  u(r,t)  =  c(r,t),  etc.  (A.2-3b) 

where  the  latter  represents  a  "hard"  boundary  and  the  former  represents 
comparatively  small  changes  in  the  continuous  medium.  The  Discontinuum 
can  be  localized,  e.g.,  particles,  bubbles,  etc.,  or  distributed,  like 
the  air/water  interface.  It  is  then  seen  that 


oM | Inhomog 


'  B:  doM  Discont. 


(A. 2-4) 


where  Qcoe.ff  is  the  coefficient  of  p(r,t)  in  the  Langevin  equation,  e.g., 
of  c(r,t)  in  (3.2),  viz.,  c~“(3^/3t^),  etc.  In  consequence  of  the 
assumed  lack  of  motion  of  the  fluid  as  a  whole,  there  is  no  doppler 
associated  with  u (_r , t ) ,  whereas  there  is  always  a  (random)  doppler  asso¬ 
ciated  with  the  Discontinuum  elements,  about  some  equilibrium  position 
of  the  ocean  wave  surface,  or  particle  displacement  in  the  volume.  (For 
a  general  treatment  of  these  important  doppler  effects,  see  [36].) 

Finally,  to  anatomize  the  various  types  of  possible  scattering  inter¬ 
action,  we  invoke  the  following  Basic  Decomposition  Principle  (BDP)  for 
the  number  of  available  dse's,  dN,  [l],  [2].  The  BDP  states  that: 


"From  the  viewpoint  of  radiative  interactions  among  the  dse's 
which  constitute  the  medium, the  argument  (dN,  cf.  (A. 2-1)),  of 
the  "counting"  functional  (/^dt^  )  =  N,  (A. 2-la)),  can  be  decom¬ 
posed  into  k  =  l,  2,  distinct  subsets  of  available  dse's  for 

radiation,  each  subset  of  which  contains  respectively  distinct 
and  independent  k-couple  radiating  elements  only."  (A. 2-5) 

This  may  be  expressed,  after  some  manipulation,  in  terms  of  (radiative) 

(k) 

mean  values,  e.g.,  <  >R,  and  fluctuations  of  the  various  densities ,  vv  ' , 
of  different  classes  (k>0)  of  dse's,  by 


>(r,t)  =  l  v'k ' (r ,t)  =  <v(r  ,t)  ^  >R  +  I  Avtk)(r,t)  ; 
k=l  R  k=l 


*>R  ’  ^ 


Av(k)  =  v(k)  _  <x;(k)>R, 


(A. 2-6) 


since  N6  =  v  ( t ) ,  r  e  A^,  when  M=S  or  V  (or  B),  of  course.  Here  (k  =  0) 
represents  the  average  number  of  dse's  (discrete  or  continuous)  engaged 
in  (re-radiation);  k=l  denotes  the  aggregate  of  independent,  uncoupled 
dse's,  while  k>2  are  the  sets  of  pairs  (k  =  2),  triples  (k=3),  etc.  of 
dse's  available  for  (multiple)  scattering.  It  follows  directly  from 
the  BDP  that  all  sets  (k>0)  are  statistically  independent.  From  this 
it  can  be  shown,  furthermore,  cf.  [21],  that  the  fluctuations  Avv  ;  obey 
(zero-mean),  k— -order  poisson  statistics,  under  very  general  conditions, 
reducing  to  gaussian  statistics  whenever  the  process  densities  are  suf¬ 
ficiently  great.  Figure  A. 2-1  gives  a  schematic  illustration  of  the 
decomposition  principle  presented  here. 

A, 2-1.  Operator  Structure  and  Radiation  Event  Statistics 

Further  insight  into  the  nature  of  the  BDP,  (A. 2-5),  and  the  canonical 
structure  of  the  inhomogeneity  operator  Q^,  (A. 2-1),  may  be  gained  from  the 
following.  It  can  be  shown  for  small  regions  dA^  in  which  dQ^  is  accord¬ 
ingly  the  scattering  operator,  that  the  BDP  can  be  anatomized  as  follows, 
cf.  Fig.  A. 2-1: 


IS)?]:  dQ  .  7  dQ<k>  !  I  a(kMkMk>:-, 

k=l  a"  k=l  0 


(A . 2-7a ) 
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the  k— -order  "quasi-particles,"  and  an  appropriate  weighting  function, 
noted  below.  Thus,  the  BDP  shows  how  a  scattering  region  is  decomposable 
into  independent  groups  of  interacting  elements.  [This  decomposition  has 
certain  analogies  with  "Cluster  Theory"  in  the  statistical  mechanics  of 
a  classical  nonideal  gas  [42],  but  in  any  case  we  are  not  concerned  here 
with  the  details  of  the  k-tuple  interactions.]  Each  k-tuple  is  entirely 
equivalent  to  a  k— -order  (classical)  "quasi-particle,"  in  the  language 
of  particle  physics  ([41],  Chapter  2).  Thus,  (A. 2-6)  can  be  interpreted 
as  a  sum  of  different  densities  of  "quasi-particles,"  which  interact  with 
the  incident  field  to  produce  corresponding  orders  of  single-  and 
multiple-scatter. 

Comparing  (A.2-7b)  and  (A. 2-1)  specifies  <£.  ,  namely, 

aM 


mi-  Q^fdt*  /  l  v(k)Y^)hfk)(t-tl,t^1,r,...)do(r,t|L',t')[  ]r, 
-oo  A,,  k-1  —  ’ 

(A. 2-8) 


t,dAM 


(k)  (k) 

where  a'  -  h'  '  is  a  time-varying  weighting  function  representing  the 

I  £ 

reradiation  process  from  the  k - order  "quasi-particles"  (or  scattering 

A 

elements*).  Again,  dQ  is  a  local  interaction  operator  (A. 2-la),  deter¬ 
mined  by  the  particular  langevin  equation  (3.2),  for  example  (A. 2-4)  et  seq. 
The  key  result  of  the  BOP  is  that  it  removes  our  having  to  specify  the 

1.  L. 

details  of  the  radiation  interactions  within  the  k— -order  "quasi- 
particle."  Instead,  we  have  a  linear  combination  of  quantities  (~y^  'v  '), 
cf.  (A. 2-8),  which  are  to  be  determined  explicitly  from  the  appropriate 
physics  of  the  scattering  situation  in  question. 

At  this  point  we  consider  further  the  concept  of  "radiation  events" 

(Z^)  (and  their  statistics),  i.e.,  the  interactions  of  the  incident  field 
with  a  quasi -parti cl e ,  followed  by  emission,  for  given  space-time 


*The  vectors  jr ,  jr]  are  defined  in  a  coordinate  system  based  in  which  is 
generally  different  from  the  coordinate  system  relating  source  (T),  the 
"ideal"  observer  in  the  field  (at  P(F[,t)),  or  the  receiver  (R),  cf.  Fig. 

A. 2-5. 
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fixed.  The  classical  description  of 


coordinates  (S=(r.,t)),  i.e., 
scattering  (Z  ,  )  is  thus  analyzed  into  a  two-component  process: 

*“~C  1  d  S  S 

Z  ,  =  l-  xZn,  and  the  classical  statistical  averages  are  <>  ,  „  = 

*-CldSS  ClaSS 

<  >  =  «  >R>5Xg*  where  <  >R  =  average  over  radiation  events  (Z^),  and 

<  >SxQ  =  average  over  spatial  and  parameter  values. 

The  importance  of  this  decomposition  is  that  it  permits  us  to  deter- 

(k) 


mine  the  overall  statistics  of  the  v' 


This  is  accomplished  by  our 


first  noting  that  the  fluctuation  Av^  in  the  number  of  k— -order 
radiation  interactions  'in  A^  is  given  by  (A. 2-6). 


Furthermore,  we  have  the  identity 


l 


(k)  (k) 

Y  v 


k=l 


h(k)  = Y(0)h(°)<v(°)>  +  £  Y(k)h(k)Av(k)  _  r(0)  +  ^  Y(k)h(k)Av(k)^ 
F  _ S  k=l  0  h  k=l  0 


k  =  0 


(A. 2-9) 


so  that  by  (A.2-7b) 


QM  =  ^A(r(0)+  I  Y^k)4k)Av(k^) 
w  M  k=l  ° 


(A. 2-9a ) 


It  can  be  readily  shown  [11]  that  the  radiation-event  statistics  (RES) 
of  the  Av^  ar( 
function  (c.f . ) 


of  the  Av^  are  zero-mean  poisson,  with  the  n— -  order  characteristic 


F(k)(H,U|ZSxa|4v(k))  -  w[/A  .  J 

(each  kM,  any  nM);  A^  n  =  •  •  *An»  (A. 2-10) 

in  which  i£=  i  -,y. . .  ,iCn;  it,  =  i  t^ , . . .  ,i  tp  and  R^  =  ,tj  ;  • . . ;  ,tp )  >  0 , 

(k,n>l)  is  the  associated  process  density.  Since  the  radiation  events  for 

different  k  are  independent,  the  complete  (conditional)  c.f.  of  the 
.  (k) 


,  ks*  1  is 


:  (°°)  . 


(k) 


k=l 


(A. 2-11) 


( k )  ( k )  ( k )  th 

The  overall  statistics  of  yv  'hr  Av  '  are  found  from  the  (n - order 

/  i.  \0  >  «  /  i.  \  M.  \  /  i.  \ 


(kl 

Fourier  transform  of)  <FV  '>7  ,  and  for  \ 

n  ZS->  1 


.  ( k ) h ( k ) ,.  (k) 


ax 


from  the  F.T. 
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{<F^  ^>,  }.  Thus,  a  central  problem  in  applying  our  canonical  development 

n  ^Sxfl  (kl 

here  is  to  determine  the  process  densities  fr  '  for  the  specific  physical 

problem  in  question.  We  note  some  examples  below  is  Sec.  A. 2-2. 

A. 2-2  Examples:  FOM  and  Classical  Scatter  Models 
We  begin  by  considering 

A.  FOM  Models: 

As  noted  in  Section  5,  FOM  models  [28]-[30]  are  quasi-phenomenological 
in  that: 

(i).  they  do  not  explicitly  include  boundary  conditions,  but  rather 
subsume  them  in  the  F  ^  and  factors  of  the  scattering 

operator,  QM>  cf.  (A. 2-9); 

(ii).  they  permit  only  k=l  (and  hence  k  =  0)  types  of  radiation  inter¬ 
actions,  namely,  the  so-called  independent  "poinf'-scattering. 

This  latter  constraint  ensures  only  "weak-scattering"  in  volumes 
(M  =  V) ,  and  correspondingly  for  surfaces  (S,B),  e.g., 

a  ;  °„  +  W H 

in  (3.9),  cf.  (5.1)  also,  with  the  added  condition  that  all 
k-order  (k^2)  interactions  in  are  omitted:  no  multiple 
scatter,  cf.  (5.8). 

Thus,  in  the  FOM  cases  we  have  the  following  relation  for  the  governing 
process  density  : 

Rlk)  "  "n1’  s  aM1)(l)j25<I* -X  e-i>'  <A-2'12> 

where  now  =  average  density  of  scatterers  in  the  illuminated  region 

(M=S,B,V).  From  (5.9),  (5.9a)  we  see  that  <v^>  =  v^=0, 

k>2. 

The  original  FOM  models  for  surface  scatter  [28] - [ 31 ]  postulated  a 

( k  1 

flat  surface  ('.  =  0),  but  employed  a  dopplerized  path  delay,  t^ ,  in  hp  ' 
which  implicitly  contained  the  effects  of  the  moving  surface.  FOM  model 
extensions,  F0M^a^ ,  explicitly  include  the  surface  elevation  ~  in  the 
doppler  delay  [36],  viz. 


A- 1 1 


-tz 


Figure  A. 2-2  Geometry  of  bistatic  surface  scatter. 


ld  ■  To +  T0  =  (ror  +  rot)/co; 

(A. 2-13) 


C  =  Jz5’ 

•v-  L 


where  r =  i  x  +  l  y  on  Ac  U  =  0)  and  iT  =  RT/  RT  ,  etc. ,  with  c  =  mean 
-  ~x  -y"  S  ~T  ~T  '-T1  0  (?) 

speed  ot  sound,  as  shown  in  Fig.  A. 2-2.  A  second  extension  ( FOM v  ;)  of 
the  original  FOM  model,  in  addition  to  the  above  doppler;ized  delay  (A. 2-13), 
explicitly  includes  the  tilted  surface  T.  of  surface  wave  heights  in  the 


"cross-section"  -y 


(1) 


e.g. 
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(A. 2-14) 


(1)  (l)~  \  (l)''  /  ~  O  \ 

n  _Yn  Jdg*  ~T  _v1r)  *  JQq*(IoT  "  loR^* 


g  ~oT  ~ oR J 


where  now  n  is  the  normal  to  the  (large-scale)  component  of  the  gravity- 
“^3 

capillary  wave  surface 


V< 


it  +  i  r  -  i 
~x^x 


'J 


(  /l 


r  v  ~z  ^ 

2 
y 


+4+? 


«Lc 


xL,x  +~y^y 


•  i  )n  }  . 
~z  z’ 


(A. 2-15) 
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For  both  extensions  r  =  a  (r)v  ,  where  again  0  is  the  mean 

0  f la )  (2) 

density  of  (independent)  scatterers  on  the  surface.  FOf'T  '  and  FOM^  ' 

thus  incorporate  explicitly  the  physical  geometry  of  the  (wave)  surface. 

The  practical  utility  of  FOM  models  is  that  they  replace  the  often 
very  difficult  problem  of  explicitly  evaluating  the  boundary  conditions 
at  the  scattering  elements,  by  introducing  a  time-variable  weighting  func¬ 
tion  yi^h!^,  or  yi^hl^,  which  also  may  include  the  extended  doppler 
or  or  Mini 

delay  (A. 2-13)  [36].  Simple  choices  of  y^  ,  h^  '  are  then  tested  against 

experience  [43]-[46].  The  limitations  of  the  FOM  models  are  threefold: 

(1) .  They  lack  explicit  physical  structures.  [This  can  be  to  a 

(2) 

considerable  extent  overcome  by  using  FOM^  ;  models,  cf.  above]; 

(2) .  They  are  quasi-phenomenological ,  namely,  they  require  a  suitable 

choice  of  h^  and  consequent  calibration  to  the  problem  at  hand;  and 

(3) .  They  (currently)  omit  multiple-scatter  effects,  or  diffraction 

f  k  1 

(e.g. ,  vv  ’  =  0,  k^  2).  [This  is  also  a  limitation  of  most  "classi¬ 
cal"  scatter  models,  as  well  [25]-[27],  [37]-[39],  [48],  but  see  [51].] 


B.  Classical  Scatter  Models: 

Classical  scatter  models  cf.  [25]-[27],  [48],  as  noted  in  Sec.  5 
earlier,  in  their  practical  formulations  for  the  most  part  neglect  multiple 
scatter  effects  (k>2),  so  that,  in  our  formalism,  cf.  Appendix  A.l, 

=  =  k^2.  In  the  important  cases  of  surface  scatter,  for 

example,  it  can  be  shown  on  comparing  the  classical  approximations,  namely 
Kirchoff  and  perturbational  boundary  evaluations  [48],  that  now 
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(A. 2-16) 


,ov(0)  i  T0V(1)  i  esf-  S  •(iT-lR). 


cf.  (5.10)  above,  with  (A. 2-15).  Equations  (5.6)-(5.7c)  again  apply  here, 
so  that 


A  a  A 


/V  A 


Kq|  surface  =  ^i-T  "-R^l  ^g*  Ij  “  JrV  '  <(7>i<(^>2  (A-2“l7) 


for  the  covariance  of  the  incoherent  scatter  from  a  typical  random  (wave) 
surface.  A  detailed  evaluation  of  (A. 2-16)  is  described  in  [50],  Part  I, 
and  extended  to  include  explicit  doppler  measures  in  Part  II.  Higher-order 
moments  are  found  similarly,  viz. 

<QlQ2Q3>surface  =£{  <Yo1^{1^02v21)yo3v3^>}  ’  etC”  (A. 2-18) 


with  the  help  of  (A. 2-16)  explicitly,  see  [50]. 
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Appendix  A.3  Bayesian  Estimation:  A  Decision-Theoretic  Formulation 

Here  we  briefly  summarize  the  needed  main  elements  of  the  decision- 
theoretic  formulation  of  signal  estimation  theory  [15],  [58],  [20]. 

For  an  estimator  we  write 

X  =Jo(-^-)  =  9E(X)>  '  (A. 3-1) 

where  8  =  (0,,...,0£j)  is  a  set  of  signal  parameters  (or  waveform  samples, 
etc.)  to  be  estimated;  X  =  (X14...,Xj)  is  the  set  of  received  data 
samples  in  which  the  estimate  is  to  be  based;  the  "decisions,"  or  esti¬ 
mates  made  are  denoted  by  Y,  =  (Y^,. . .  ,Y£j),  and  o  in  Y0  indicates  the  a 
priori  probability  (density)  governing  the  parameters  j3,  e.g.,  o  =  o(9), 
or  S_  in  o(Sj-  [We  note  that  Ya  is  an  estimator  for  all  permitted  X, 
while  for  a  particular  set  X  =  V  ,  Y0  becomes  an  estimate.]  The  decision 
rule  6  is  here  5(Y|X)  =  £(Y-Ya(0 1 X.) ) ,  which  is  an  (M-dimensional )  delta 
function.  When  Y0  (=g£)  is  given,  i.e.,  the  estimating  receiver  gE(X)  is 
specified,  the  average  error  (or  ri sk)  is  determined  from 

R(o  ,6 ) c  nv.  0  =  /  o(S  or  0  )dS  (or  de)  •  {  F  (X  |  5(9.)  )C(S  or  0  , Y  )dX 

‘  (A. 3-2) 

where  C  is  an  appropriately  chosen  cost  function,  naturally  proportional 
to  the  measure  of  error  selected.  The  estimator  Y_  is  a  point  estimate, 

—o  i - 

and  is  the  conditional  pdf  of  X,  given  S(0). 

Another  useful  type  of  estimation  procedure  is  interval  estimation, 
defined  by  the  probability  P  that  a  particular  point  estimate,  Y  (for  given 
X)  falls  within  (1±\)  100"  of  the  true  value  of  the  quantity  [S.  =  waveform, 
or  0  =  parameters  in  S(o)]  being  estimated,  viz.: 

P{  (1-X  )[S  or  X']  ^  Y  (£  °r  0  |X)  <  (1+A  )[S  or  6]},  0  -  A  <  1 

(A. 3-3) 

where  »  is  a  measure  of  the  prechosen  confidence  interval  selected,  -AS.  to 
+>S  for  example.  In  the  case  of  a  single  parameter  0,  for  instance,  we  can 
wri  te 
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(l+x)e 

P{(1-A)0  <  Y  (0|X)  <  ( ?.+X) 9}  =  /  W  (Y  (G))dY 

(1-X)0  1  ° 


(A.3-4) 


where  W^(Ya(0)  =  p ( y | 0 )  is  the  conditional  pdf  of  the  estimator  Y,  con¬ 
ditional  on  0,  formed  here  from  p ( Y 1 6 )  =  /  Fn(X.|S.(9) 6(yq |X_)dX ,  with 
obvious  formal  extensions  to  the  multidimensional  cases,  cf.  [163,  [58], 
[20].  The  unconditional  pdfs  of  the  estimators  themselves  are  obtained 
from  p( y)  =  <p(Y|S.or  6)>s  Qr  0. 

We  emphasize  that  YQ~is  appoint  estimator,  embodying  the  specific 
structure,  g£,  (A. 3-1),  of  the  receiver  performing  the  estimation.  On 
the  other  hand,  interval  estimators,  as  expressed  by  P,  (A. 3-3),  (A.3-4), 
yield  a  probability  which  is  a  measure  of  the  efficiency  of  the  point 
estimator  for  any  particular  application  (i.e.,  choice  of  Xj .  The 
average  error  (or  risk)  (A. 3-2)  measures  the  expected  cost  or  average 
error  in  using  ja,  considered  over  all  possible  {X}  received. 

For  optimum,  or  Bayes  estimation  we  seek  estimators  Y*  which  minimize 
the  average  error  or  risk  R(o,6),  (A.  3-2).  The  general  form  of  the  re¬ 
sulting  depends,  of  course,  on  the  choice  of  "cost11  or  error  function 
C(S  or  0,Ya).  For  example,  for  the  quadratic  cost  function  (QCF) 


C(u,Yj  =  C  |0-Y  I2  =  C  1  (0~ 
ou-o'  0  c  ,  m 


m=l 


■Y  ~)‘ 
am' 


(A. 3-5) 


the  associated  optimum  estimator  is  found  to  be  (cf.  Chapter  3,  [58],  Chapter 
21,  [20])  the  set  of  equations 


io(9|X)QCF  =  L  °'[£)FnU  !S(0))dO, 


(A. 3-6) 


wi  th  o(  9) -►  o(5>) ,  dO-*dS  in  the  case  of  estimating  signal  waveforms. 

Note  that  Y*  is  generally  a  nonlinear  operator  on  the  received  data,  X. 

Another  cost  function  of  considerable  interest  is  the  simple  (or 
rectangular)  cost  function  (SCF)  given  by 


M 

1 

in=  1 


[A  (  Y 0"  )  ] 

L  m  '  m  m' J 


(A. 3-7) 
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with  appropriate  choices  of  CQ,  A-  to  ensure  meaningful  results  (e.g.,  posi¬ 
tive  errors,  etc.)*  Minimization  of  average  risk  here  leads  directly  to 
(cf.  Sec.  21.2-1  of  [20])  the  following  relations  determining  Y*| 


90' 


log  {oCejW-iQCj0-)} 


0-'=0^=Y? 

m  m  m 


=  0,  all  m  *  1 , . . . ,M,  I  (A. 3-8) 


-here  WjUIy-6;,)  =  <F„<il  S(Y;,,e' >», ,  /  afy^e,.) ,  end 


Y~ 

—  -  ~m 


°(Yfts0fi|)  =  <a  (y^  ,9J  )>fl  , ,  withQ.1  =  all  0  except  j9~.  The  condition  (A. 3-8) 


m 


determining  y*  =  (y*)  is  precisely  that  determining  the  unconditional 


maximum  likelihood  estimates  (UMLE's)  of  9  =  (9^> ,  viz.  Y*  here. 
-  *  m  m 


[If 


o(0ft)  is  omitted  in  (A. 3-8)  one  has  the  corresponding  conditional  maximum 

likelihood  estimates  (CMLE's)  of  the  {9~}.  The  extension  of  (A. 3-8)  to 

waveform  estimation  is  formally  immediate:  one  replaces  9-  by  S-,  cr(0ft) 

J  r  m  m  m 

by  o(Sft),  W^Xj0-)  by  Wj(XjS~),  etc.,  with  m  =  1,...,J  now.] 


In 

I. 


II. 


III. 


summary,  we  remark  that: 

The  maximum  (conditional)  likelihood  estimator  (CMLE)  maximizes 
the  probability  of  a  correct  decision,  without  regard  to 
incorrect  decisions  and  their  costs; 

The  maximum  unconditional  likelihood  estimator  (UMLE)  maximizes 
the  probability  of  a  correct  estimate,  when  all  possible  (signal 
or  parameter)  values  are  taken  into  account,  again  without 
specific  regard  for  incorrect  decisions  and  their  errors; 

The  optimum  quadratic  estimator  Y^ ,  (A. 3-6),  is  an  uncon¬ 
ditional  least  mean-square  error  (LMSE)  estimator,  which 
accounts  on  the  average  for  incorrect  decisions  and  errors. 


N\ 

i 


t 


) 


(A. 3-9) 
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